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Abstract 

A way of construction of semiclassical wave function (SWF) based on the Maslov - 
Fedoriuk approach is proposed which appears to be appropriate also for systems with 
chaotic classical limits. Some classical constructions called skeletons are considered. The 
skeletons are generalizations of Arnolds' tori able to gather chaotic dynamics. SWF's are 
continued by caustic singularities in the configuration space rather then in the phase space 
using complex time method. The skeleton formulation provides us with a new algorithm 
for the semiclassical approximation method which is applied to construct SWF's as well as 
to calculate energy spectra for the circular and rectangular billiards as well as to construct 
the simplest SWF's and the respective spectrum for the Bunimovich stadium. The scar 
phenomena are considered and a possibility of their description by the skeleton method 
is discussed. 

PACS number(s): 03.65.-w, 03.65.Sq, 02.30. Jr, 02.30.Lt, 02.30.Mv 

Key Words: Schrodinger equation, semiclassical expansion, Lagrange manifolds, classical 
trajectories, chaotic dynamics, quantum chaos, scars 



1 Introduction 



The semiclassical approximation is widely used and by this it is the well known method of 
approximation in quantum physics. There are two basic formulations of the method - the one 
based on the wave function formulation of quantum mechanics [1, 6] and the other - on the 
Feynman paths integral [2, 3]. While both the formulations of quantum mechanics are known 
to be equivalent the general wave function formulation of the semiclassical approximation [6] 
is considered to be not applicable in higher dimensional quantum problems which classical 
limits are chaotic. 

In fact it is a common convince that the only way to formulate the method in the last 
cases is the Gutzwiller approach based on the Feynman paths integral [4]. This convince 
follows from also a common believe that the wave function formalism can be applied only in 
the cases when the classical limit of the quantum problem is the integrable one i.e. if the 
classical motion is set on the Arnold tori [5, 7] on which the semiclassical wave functions 
(SWF) are constructed. As the main argument for such a believe the KAM theorem [5, 7] is 
invoked which claims that the Arnold tori structure of the classical phase space disappears 
if classical systems become non-integrable. It is argued that because of that the only phase 
space finite structure which still survive in chaotic motion of the classical system are periodic 
orbits which therefore provide us with a skeleton on which the Gutzwiller formula is built. 

However one can criticise the point of view that the existence of the Arnold tori structure 
of the phase space is necessary for a possibility to construct the semiclassical wave function, 
i.e. that Arnold's tori provide us with a unique support for such a construction. Such a 
definite conclusion does certainly not follow from the original local approach of Maslov and 
Fedoriuk to construct SWF's [6]. 

On the other hand while the results provided by the Gutzwiller method are very rich and 
appreciated, particularly when energy spectra of the chaotic systems are considered (see for 
example [8]), the method itself does not allow us for constructing and discussing properties 
of the wave functions involved in the problems considered. In fact the wave functions are 
found in such cases by different methods mostly numerically. It is just due to such numerical 
calculations of the wave functions [9, 10] that a phenomenon of scars has been discovered [14] 
which existence in the wave function patterns is still waiting for its full explanation [16, 19, 15] 
(see also other papers cited in [16]). 

Billiards while a non-analytic motion area are well known however as examples of the non- 
integrable two dimensional systems except the known cases of the elliptical and rectangular 
billiards. They are widely considered as a simple field of experimental [17, 18] as well as 
theoretical [20, 21, 22, 23] (and papers cited there) and computational investigations [9, 
11, 12, 13] allowing to apply many different methods (see Sarnak's lecture [24] and [25] of 
the same author for an extensive review of the respective theoretical methods covering also 
billiards manifolds). 

In this paper we are going to develop the SWF formalism which can be applied at least 
in principle to non-integrable cases of the two dimensional motions in billiards and which can 
be easily extended to higher dimensions. 

Essentially our approach is initially very close to the one of Maslov and Fedoriuk [6] . The 
main difference between them is in a treatment of crossing the singular points of the SWF's 
set on caustics. Namely, instead of making the canonical phase space variable transformations 
accompanied by the Fourier transformations of the SWF's to move through the caustic points 
we apply the analytical continuation on the complex time plane to both the SWF's and the 
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classical trajectories. This greatly simplifies the corresponding procedure in comparison with 
the Maslov and Fedoriuk treatment. It is the exceptional role played by the time variable in 
the semiclassical limit of the Schrodinger equation which permits us for such simplification. 
Because of this the SWF can be considered as depending effectively on the time variable only 
while the remaining variables plays the role of spectators. 

Therefore the SWF's are first constructed locally on so called bundles of rays to satisfy 
vanishing boundary conditions. Next they are matching to get a global semiclassical solution. 
This is done however with the help of an earlier constructed set of reversible in time closed 
connected ray bundles called bundle skeleton which play a role of Arnold's tori except that 
a number of ray bundles in the skeleton can be infinite. 

This is just the notion of the ray bundles which allows us to catch a possible chaotic motion 
in the billiards not to resign from considerations of a set of trajectories on which SWF's can 
be defined while the bundle skeleton idea allows us to close the matching procedure of SWF's 
constructed locally. 

The paper is organized as follows. 

In the next section the Maslov - Fedoriuk method of the semiclassical wave function 
construction is reminded and discussed. 

In sec. 3 the semiclassical wave function is considered as the classical objects which time 
evolution is described by the classical equations of motion. 

In sec. 4 the construction of global SWF's in billiards is given. 

In sec. 5 the circular billiards is considered to demonstrate how the method works in the 
case of the presence of caustic. 

In sec. 6 the rectangular billiards is considered as the case deprived of a caustic. 

In sec. 7 our method is applied to the Bunimovich stadium to show its usefulness in 
describing the so called bounces ball modes. 

In sec. 8 we discuss a possibility to describe by the skeleton idea the scar phenomenon 
considering such a scar formed around the horizontal periodic orbit in the Bunimovich sta- 
dium. 

In sec. 9 the results of the paper are summarized and some limitations of the method are 
discussed. 

There are four appendixes attached to the paper which justify the main assumptions used 
in the construction of the global SWF's on skeletons. 

2 Semiclassical wave function expansion for n-D stationary 
Schrodinger equation 

Consider the n-dimensional stationary Schrodinger equation: 

A^(r) + A2^(£;-y(r))*(r)=0 (1) 

with a potential V(r), r G Rn confining a point particle with a mass m and containing a 
formal dimensionless parameter A > 0. For a convenience we shall put further h = 1 and 
m = 1. The Schrodinger equation is recovered by putting A = 1 in (1). 

We would like to construct a solution to Eq.(l) using the idea of Maslov et al [6] and 
considering the wave function ^'(r) as defined on families of classical trajectories a dynamic 
of which is given by the classical Hamiltonian H = + V{r) and which carry an energy 
Eq all. 
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Such families are constructed locally in the following way. In i?„ we choose a n—l-D hyper- 
surface Sn-i on which the initial momenta p(ro), tq G Sn-i) will be defined so that the pair 
(ro,p(ro)), ro G will serve as initial data for the trajectory r{t) = f (ro, p(ro); t). Then 

the transformation: r ^ (t, ro) ^ (i, si, Sn-i), r = (xi, x„), tq = (xo,i(si, Sn-i), 
2^o,n(si, Sn-i))) ((^ii ■5n-i) parameterize the hypersurface is one-to-one up to a 

caustic surface C„_i on which the Jacobean (f (ro, p(ro); = f (i, si, 



dt ' dsj 



(2) 



vanishes. 

A n-dimensional domain A„ of 2n-dimensional phase space i?2n made in this way by the 
hypersurface and trajectories emerging from it is known as the Lagrange manifold [5]. 

Therefore in the variables t, si, the new wave function x(i, si, satisfies 
the following relation with the previous one: 

|x(t,si,-,s„_i)p = |'I'(f(t,si, ...,s„_i)p|J(t,si, ...,s„_i)| (3) 

The particle momentum p on the trajectories v{t) = f (t, si, s„_i) satisfies of course 
the equation: 

df{t,Si,...,Sn-l) 1^,. 

=p(f(t,si,...,s„_i)) (4) 

defining also the Jacobean evolution. Namely: 

d dfi{t,Sl,...,Sn-l) _ ■A dpi dfj{t,Sl, ...,Sn-l) 

dt dt ~^,dx-j dt 

■1=1 J 



d dfiit, Si, Sn-l) _ ^ dpi dfj{t, Si, Sn-l) 



- V 

dt dsi r— ^ dxj dsi 



so that 



Z = l,...,n-1 (5) 

^'^''^'^^gt '^""^^ ^ J(t,si,...,s„-i)Vp(f(t,si,...,s„_i)) (6) 
The above equation is just the Liouville theorem with the solution: 

J(t,.i,. = J(.i,...,5„_i)e/o^P(f(*'^^--^"-))'^*' (7) 

where J(si, is the value of the Jacobean on the hypersurface 

It is well known from the classical Hamiltonian mechanics [5] that the action integral: 

5(r,ro)= rp{r')dr' (8) 

Jro 

taken on the Lagrange manifold A„ is a point function of r and ro- Therefore taking ro 
as a definite fixed point of the hypersurface Sn-i and denoting by S{r) the action function 
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corresponding to this case we can complete a definition of the wave function x{t, si, s„_i) 
by the following equation: 

*<^(f (t, Si, Sn-l)) = J--^ {t, Si, S„_i))e<^^^^(f(*'^l--^"-l)x'^(t, Si, Sn-l) 

(9) 

where cr = it is a signature of *'^(r). 

Therefore the quantities involved in the above definitions satisfy the following equations: 

p(r) = V5(r) 
^p2(r) + y(r)-i?o = 

A( J-^X'^(r)) + a2iXJ-^r)Vx''{r) ■ p(r) + 2X'^{E - £;o) J'^ (r)x''(r) = 

r = f(i,si,...,s„_i) (10) 

By the variables t, si, Sn-i the third of the last equations can be rewritten in the 
following form: 

"^'^ dt + 

J5A(j-5x'^(t,si,..,s„_i,A)) +A'(E-£;o)x'^(t,si,...,Sn-i,A) = (11) 

where a dependence of x'^(t, si, A) on A was shown explicitly. 

The Eq.(ll) describes the time evolution of x'^(^) "^i, -Sn-i, A) along trajectories starting 
on the hypersurface if its "initial" values on this surface, i.e. x^(0) ^i) •■•j ^n-i) -^) = 

X^isi, Sn-i, A) are given. 

We are going to consider the equation (11) in the semiclassical limit A +oo looking for 
its solutions in the form of the following asymptotic series: 

E — Eq = ^ E^X ^ ^ 

k>l 

X''{t,Si,...,Sn-l,X) = ^XUi^Si,...,Sn-l)X~'' 

k>0 

X"(S1, Sn-l, A) = ^ Xkisi, Sn-l)X-'' (12) 
fe>0 

Putting A = 1 in (1), (9) and (12) we get approximate semiclassical solutions to the energy 
eigenvalue problem of the Schrodinger equation. 

It is to be noticed that for the selfconsistency reasons the semiclassical series for the 
energy parameter in (12) starts from the second power of A^"*^, i.e. this ensures the proper 
hierarchy of steps in the algorithm of semiclassical calculations by which the higher order 
terms of the series in (12) are determined by the lower order ones. 

It should be noticed also that despite the fact that Eq enters the classical equation of 
motion (10) it is still quantum, i.e. its value depends on h which is considered to have the 
definite numerical value, i.e. h is not a parameter. In particular the series (12) represent 
the inverse power hierarchy in the formal parameter A, i.e. not in powers of h, between 
subsequent terms. 
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Moreover Eq if quantized can depend on A. However, whatever is this dependence the 
semiclassical series of the difference E — Eq must be given by (12). 

Needless to say the introducing A makes a treatment of the Schrodinger equation equiv- 
alent of course to considering it in the limit 7i — )■ 0, i.e. scmiclassically, clearly however 
separating the role of /i as a parameter from its role defining the microscale of quantum 
phenomena. 

Substituting (12) into (11) we get: 

dXo{t,Si, ...,Sn-l) ^ 

dt 

dXk+l{t,Sl,...,Sn-l) ^ 

at 

' k \ 

A [j-hUt, Sl, Sn-l)) + 2 ^ Ek-l+lXl{t, Si, Sn-l) 
V 1=0 / 

k = 0,1,2,..., (13) 



at 



with the obvious solutions: 



Xo{t, Si, Sn-l) = Xoisi, Sn-l) 
Xk+l{t,Sl,...,Sn-l) = X^+i(si,...,S„_i) + 

[j-hUf, Si, .:, Sn-l)) + 2 ^ Ek-l+lXlit', Si, Sn-l)) dt' 



1=0 

k = 0,1,2,..., (14) 



3 Semiclassical wave function expansion elements as classical 
objects 

The recurrent system of equations (13) can be considered also from the classical point 
of view as defining the time evolutions of Xki^, si, Sn-i), k = 0,1, 2, along the classical 
trajectories r(i) = i(t, si, Sn-i). Namely, for each given initial point ro(si, s„_i) and 
a given trajectory emerging from it let d denotes a distance measured along the trajectory 
from the point ro(si, s„_i) to the point {x\,...,Xn) lying on this trajectory. The set 
{d,si, ...,Sn-i) can be used as the new coordinates instead of (xi, x„). Their mutual 
relations are given by: 



d = d{t,Si, ...,Sn-l) 



d{{t',Si,...,Sn-l) 



dt' 



dt' 



(15) 



and 



r = {xi,...,Xn) = f(t(rf, Si,..., Sn-l), Si,..., Sn-l) = g(c/. Si , S„-l) 



(16) 



where t = t{d, si, Sn-i) is the solution of d = d{t, si, Sn-i) with respect to t. 

Considered on the Lagrange manifold A„ the system is then governed by the Hamiltonian: 



Ha = 



1 



Pd 



ds{d,si,...,s„-i) 
dd 



+ V{d,Si,...,Sn-l) 



(17) 
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where V{d, si, Sn-i) = V{i{t{d, si, Sn-i), si, Sn-i)) and is the momentum conju- 
gated with d. 

The corresponding classical equation of motion for x'^(t, si, A) on the manifold 
Aji is therefore: 

= {X ,i^Aj + = (18) 

since x'^{t, si, Sn-i, A) being independent of d commutes with the hamiltonian Ba„- 
Therefore, from (12) and (13) we get on the Lagrange manifold A„: 

dXoit,Si, ...,Sn-l) ^ Q 



dxUlit,Sl,...,Sn-l) ^ ^jiA(j-^X^(t,.i,...,.„_i)) +^^E^.-mXr(t,^l,-,«n-l) 

fc = 0,l,2,..., (19) 

Of course (14) provides the solution to the equations (19). 

To conclude Xo(^' ■Sn-i) is constant on the Lagrange manifold according to (14) 
being obviously independent also of d. 

It is to be noticed that defining on another Lagrange manifold by choosing a dif- 

ferent set of trajectories we obtain also another SWF defined on this new Lagrange manifold, 
i.e. different from the one defined on A„. This is why the corresponding Xo •■•) ^n-i) 
cannot be considered as the global integral of motion, i.e. in the whole phace space i?2n of 
the Hamiltonian H - each particular SWF is defined only on a particular Lagrange manifold 
A„ corresponding to it. 

For future applications it is worth to note that if (11) is obviously invariant on a reparametri- 
zation of the hypersurface it is also invariant on the following change of variables: 

t t(si, ...,Sn-l) - t 
Sk hk{si,...,Sn-l) 

k = l,...,n-l (20) 
if it is accompanied simultaneously by the transformations: 

X'^{t,Si,...,Sn-l,\) 

_ 1 

X~°'(t,si,...,s„_i,A) = Ji 2(/if^(si,...,Sn_i),...,^~ii(si,...,S„_i)) X 

X'^(r(si, Sn-l) - t, hi^{si, Sn-l), K^^iisi, S„_i), A) (21) 

where Ji(si, s„_i) is the Jacobean of the transformation (20). 



4 Semiclassical wave functions in the billiards 

Let us apply the above formalism to construct continuous semiclassical wave functions 
inside billiards B (Fig.l) vanishing on its boundary. Such a construction will be done in 
several steps the first one consisting of some geometrical preliminaries describing a notion of 
a skeleton, i.e. the closed set of families of trajectories which forms a base on which SWF's 
are constructed. 
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Figure 1: An arbitrary billiards 



4.1 Ray bundles and bundle skeletons 

We shall assume that the billiards is classical, i.e. its boundary dB is a closed curve 
independent of A and given by r = ro(s) = [xo(s), yo(s)] where s is a distance of a boundary 
point ro(s) measured clockwise along dB from some other point A of dB chosen arbitrary, i.e. 
s{A) = 0. Both xq{s) and yo(s) are continues. The curve however consists of a finite number 
q, q > 1 oi smooth arcs Ai,...,Aq with respective length Li,...,Lq, so that the derivatives 
Xq{s) and ?/o(s) are discontinuous in a finite number of points on the segment < s < L 
where L = Li + • • • + Lg is the total length of dB. Both xo(s) and yo{s) are of course periodic 
with the period equal to L. We shall identify the point A with the point beginning the arc 
Ai. 

Next we define a bundle of rays as a family of trajectories in the following way. 

Let Ai~{u, I), Li + • • • + Lk-i < u < Li + ■ ■ ■ + L^, < I < L^, be an open connected 
piece of the arc Ak beginning aX s = u and having a length I. 

Let further Tk{s,t]u,l), u < s < u + I, < t, be a family of trajectories given by 
angles 7fc(s;n, < 7fc(s;u, /) < 27r, at which the trajectories escape from Ak{u,l). The 
angles 7(5; n, /) are smooth functions of s and are measured with respect to the x-axis while 
the tangential vectors t(s) = [ "^^^g^^ ; ^^di^^ ] ~ [cos /3(s), sin /3(s)] are inclined to the x-axis 
at angles /3(s) (Fig.l). The latter angle can be discontinuous at the points where Xq(s) 
and y'Q{s) are discontinuous. Then the angle ak{s;u,l) = 7^(5; u,/) — (3k{s) is made by 
the classical ball momentum p(s;u, /) on the trajectory with the tangent vector t(s), i.e. 
p(s; u, I) ■ t(s) = pcos ak{s; u, I). It is assumed that < ak{s; u, I) < tt. 

The classical time evolution of the family rk{s,t;u,l), u < s < u + I, is therefore the 
following 

rkis,t;u,l) = ro{s) + p{s;u,l)t, ro{s) £ Ak (22) 
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where p{s;u,l) = [pcos^k{s',u,l),psmjkis',u,l)] satisfies the classical equations of motions 
(10), i.e. p'^(s;u,l) = 2Eq (again we put m = 1 for the billiard ball mass). 

The trajectories (22) define of course the change of variables {x,y) — >■ {t, s), in vicinity of 
Ak{u, I), i.e. X = fk{t, s; u, I), y = gk{t, s; u, I) with the Jacobean: 

Jk{t, s; u, I) = p^j'ki^; u, l)t — p \t{s)\ sin 0:^(5; u, I) = p'^^^i^; u, l)t — psin Qfc(s; u, I) 

(23) 

since |t(s)| = 1. 

The family of trajectories defined in the above way will be called a bundle of rays 
emerging from the segment Aj^{u, I) of A}^ and will be denoted by Bk{u, I) while the trajectory 
itself will be called rays. 

Since each ray of the bundle B^iu, I) after some time rfc(s; u, I), ro(s) € Ak{u, I), (different 
for different rays) achieves another point of the boundary dB it means that the bundle Bk{u,l) 
maps the segment Ak{u,l) into another piece BAk(u,l) of the boundary dB. In general it 
is assumed that this map of Ak{u,l) into BAk{u,l) provided by the transformation (22) is 
one-to-one except the caustic points of BAk{u,l) in which JK{Tk{s;u,l),hk{s;u,l);u,l) = 
0, ro(/iA;(s; n, /)) G BAk{u,l). Here hk{s;u,l), ro(s) G Ak{u,l), realizes explicitly this map. 
The caustic points will be assumed to be isolated singular points of the map Ak{u,l) — > 
BAk{u,l). Typically they need a special care if they are. 

By DBk{u,l) will be denoted a domain of the billiards B covered by rays of the bundle 
Bk{u,l) emerging from ^^(tt,/) and ending at BAi;[u,l). The domain DBk{u,l) is locally a 
Lagrangian manifold on which each loop integral § p ■ dr vanishes. 

Assume further the boundary dB to be a mirror-like, i.e. reflecting the incoming rays 
according to the reflection principle of the geometrical optics and let Aii.i(u',l') be a piece of 
another arc A/;/ of dB such that Ak'{u', I') D BAk{u, /) 7^ on which another bundle of rays 
Bk'{u', I') = {rjfc/(s, t; u', V) : u' < s < u' + I', < i} is defined. 

If on the segment Ay{u',l') n BAk{u,l) the ray bundle By{u\l') coincides with the 
reflected one we call the ray bundle Bk'(u',l') a reflection of the bundle Bk{u,l) on the 
segment mentioned. 

The reflection operation over the bundle Bk{u,l) will be denoted by 11 so that IlBk{u,l) 
denotes the set of all rays arising by the reflection of the bundle Bk{u, I) of dB. 

Consider now a family of the disjoint ray bundles B = U Bk{u, I), B^iu, l)r\Bki(u', I') = 
ifBk{u,l)^Bk'{u',l'). 

The family B will be called closed under reflection 11 on the boundary dB if the following 
two conditions are satisfied for each bundle Bk{u,l) G B: 

n 

UBk{u,l) = [j Bj{uj,lj)r\UBk{u,l), Bj{uj,lj) G Ai., j = l,...,n 

m 

Bk{u, = U ^Bi{ui, k) n Bk{u, I), Bi{ui, k) G Aj^, i = 1, m 
1=1 

HB = B (24) 

A closed family B' = [J Bk/{u', I') is embedded into a closed family B = (J B^i^u, I) if each 
ray bundle of B' is a subset of some ray bundle of B and each bundle of B contains some 
bundle of B'. 
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A closed bundle family is ealled eonnected if it cannot be decomposed into a sum of 
another two disjoint closed bundle families. 

A closed connected bundle family will be called a Lagrange bundle skeleton or simply 
a skeleton if it cannot be embedded into another closed connected bundle family. 

Prom now on all considered bundle families will be assumed to be skeletons. 

A basic property of each ray belonging to a skeleton B is that by its time evolution and 
bounces on the billiards boundary it will never leave B. 

Another basic property of a skeleton B is its completeness, i.e. none bundle can be added 
to or removed from the skeleton not destroying it. 

Let us now reverse in time all trajectories belonging to B. This operation leads us again 
to some skeleton B"^ which will be called associated with B. 

Bundles of B"^ are obtained simply from the corresponding bundles of B. Namely with 
each bundle Bk{u,l) of B let us associate a bundle B^{u,l) which trajectories satisfy the 
following condition: 

7j^(s;u,Z) =7r + 2^(s)-7fe(s;u,/), u<s<u + l (25) 

i.e. these trajectories are just the reflections on dB of the time reversed trajectories defined 
by 7fc(s; u, I) and belonging to Bk{u, I). 

The skeleton B"^ is organized by all bundles B^{u,l) obtained in the above way. 

4.2 SWF's defined on ray bundles of a bundle skeleton 

Consider a bundle skeleton B. On each of its ray bundle Bk{u, I) we can now define the 
following pair of SWF's s; u, I; A), a = ±: 

n{t,s;u,l-X) = j;^(t,s;u,Oe'^^^(^'*+^/o''-«"'=(^;-'0<^^');^-(t,s;u,Z;A) (26) 

where = 2Eq and X%{t-, s; «, Z; A), a = it, are given by (12) and (14). 

Exactly in the same way we can define a pair *^.j.(t, s; u, Z; A), a = it, of SWF's on the 
corresponding associated bundle B^{u,l): 

n;kit,s;u,l;X) = J2l,it,s;u,l)e'''^ip"'-Pfo-°'->'(-'-'l^^^^^ (27) 

It will be also convenient for further considerations to substitute the time variable t by 
the distance variable d = pt and consequently to give the trajectories (22), the Jacobean (23) 
and the wave function (26) the following forms: 

rk{d, s; u, I) = ro(s) + d(s; u, I) 
d(s; u, I) = pt = [dcos7fc(s; u, l),dsm 7fe(s; u, I)] 

ro{s) e Ak{u,l) (28) 

and 

Jk{d,s;u,l) = ^Jk(t,s;u,l) = ^Jjii^^hlld- smak{s;u,l) 

rois) e Ak{u,l) (29) 
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and 



ro{s)eAk{u,l) (30) 

where Xki^: s;u,l; X) = pXfc(|) •s; ^; A), a = ±, and k = pX is the wave number of the 
billiards ball. 

Nevertheless, for simplicity of notations, the bar over Xki^, s;u,l; X) will be dropped in 
our further considerations. 

By the variable d the solutions (14) can be rewritten in the form: 

Xk,o{d,s;u,l) = Xkfi{s;u,l) 
Xk,j+iid, s; u, I) = Xkj+iis; u, I) + 

iDk{a,s;u,l)xk,j{a,s;u,l) + 2 ^ Ej_„+iXfe,^(a, s; u, Z) | da 



at 
2^ 



m=0 

J =0,1, 2,..., (31) 



where Dk{d, s; u, I) = (d, s; u, l)-Ak{d, s; u, l)-J^ ^ {d, s; u, I) and Ak{d, s; u, I) is the Laplacean 
expressed by the variables d and s corresponding to the Bk{u, Z)-bundle. 

^%id, ,s; u, I; A), a = zt, arc defined initially in the domain Di^{u, I) of the billiards which 
boundary dDk{u, I) contains of course Ak{u, I). A pattern of the remaining part of dDk{u, I) 
depends on the corresponding caustic K^{u, I) = {{fk{Ki^{s; u, I), s; u, l),gk{K]^{s; u, l),s; u, I)) : 
Jk{Kk{s;u,l),s;u,l) = 0, ro(s) G Ak{u,l)}. It follows from (23) that 

Kk{s;u,l) = ^^^^^^ ^\ ro{s)eAk{u,l) (32) 
ds 

and Kk{s;u,l) can escape to infinity if ^^^^gHiO ^ q. 

If the set Kk{u,l) fl S is not empty then it belongs to dDk{u,l). The remaining rays of 
Bk{u, I) which are not tangent to K^iu, l)r\B end at BA]^{u, I) and complete in this way the 
domain Dk{u, I). 

An important property of the representation (26) is its uniqueness, i.e. for two different 
bundles defined on the segment Ak{u,l) this representation provides us with two different 
pairs of ^%{d, s; u, I; A), c = =b. This conclusion follows from the fact that for A sufficiently 
large the SWF's arc determined only by exponentials and the latter are different at the same 
points {x, y) for different bundles. 



4.3 Propagation of ^^(rf, s; I] X) through a skeleton 

^%{d, s; u, I; A), o" = it, being defined in Dk{u, I) G Bk{u, I) can be continued on the whole 
skeleton by accepting two rules the solutions have to satisfy. 

A need for the first rule to be formulated arises when the solutions propagate by the 
caustic where they are singular. These are branch points singularities on the d-plane of the 
type (c? — dc{s))~^ where Jk{dc{s), s;u,l) = 0. These points have to be avoided somehow 
when continuing the solutions along the rays of Bk{u,l). Relying on the results of App.B 
we accept the rule of avoiding these points clockwise by ^^{d, s; u, I; A) and above them and 
anticlockwise by {d, s;u,l; X) and below them. 
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The second rule has to be given to describe bouncings of "^%{d, s; u, I; A), cr = it, by the 
billiards boundary. 

Let Bk'{u',l') be a reflection of the bundle Bk{u,l) i.e. Bki{u',l') C IlBk{u,l). Continua- 
tions of s; u, I; A), a = ±, onto Bki{u' , V) is given by the following formula: 

^''^'^\d,s-u',l'-X) = -r],J~^d,s;u'j') x 

^ai\p(d-J^, cos af^,{s';u',l')ds'+Sk{u,l))^<T,cont^^ ^. ^/ ^-j 

Sk{u,l) = D{hf. {s;u,l);u,l) + / cos ak' {s' ; u' ,l')ds' — / cos ak{s';u,l)ds' 

Ju' Ju 



Xk {0,s;u,l; A) 



dhk{s';u,l) 



ds' 



Xk{D{hk^{s\ u, I); u, I), /ife ^(s; u, l)\u, I; A) 

s'=h^ {s;u,l) 

ro{s) = ro{h^\s; u, I)) + n{h^\s; u, l);u, I) (33) 



The factor rja- in (33) depends on whether there was a caustic point on the way of contin- 
uations of (rf, s; u, I; A) in which case r/o- = ai but r]a = 1 in the opposite case. 

D{s;u,l) = \Y){s;u,l)\ denotes the distance between the points ro(%(s;u, Z)) and ro(s) 
of the boundary dB. 

It is to be noticed that Sk{u,l) is s-independent (see App.C). 

x'^''^* (d, s; u' , I' ; X) is defined by (31) with the initial values given by the corresponding 

semiclassical scries of x^™" (0, s; ^i', A) as given by (33). Due to the chain property of 
the formulae (31) such a definition of Xk^°^*{d, s; u' , I'; A) ensures that Xki^^ s; u, I; A) can be 
considered as a continuous function of d on the skeleton as d changes along rays arising by 
subsequent reflections of the one starting from the bundle Bk{u, I). 

Another property of the definition (33) is that close to the arc A^./ there is a domain 
in which both the SWF's '^I'^'^^d' , s' ; u' , I' ; X) and "^^d, s;u,l; X) are defined so that the 
superposition: 

^I'^'^Pix, y, A) ^ M/^'^°"*(d', s'; u' , I'; A) + nid, s; u, I; A) 
X = xq{s) + dcos ak{s; u, I) = xo{s') + d' cos ak'{s'; u , I') 
y = yo{s) + d sin ak{s;u, I) = yo{s') + d' sin ak'{s';u' , I') (34) 

can be done. 

It follows easily from (33) that on the arc A^r, i.e. for d' = this superposition vanishes. 



4.4 SWF's vanishing on the bilUards boundary 

Another obvious property of (d, s; u, I; A), a = it, is that they cannot vanish on Ak{u, I) 
unless Xki^j s;u,l; X) vanish there identically. Therefore a wave function ^'^^''^{x,y;u,l; X) 
vanishing on j4fe(n, I) should be represented in the semiclassical limit by a linear combination 
of at least two SWF's of the form (26). The superposition (34) is an important example of 
such a combination. It is shown however in App.A on a more general level of considerations 
that the proper linear combinations have to be the following: 

^l^'^ix, y; u, I; A) = %{di, si;u, I; A) + ^^^(^2, S2; u, I; A) = 

jPidi, si; u, /)e''^*^(''i+r ^°«"'=(^''"'')'^^');^-(di, si; /; A) + 

^4(^2, S2; U, l)e-^<''^-S:" cosa,(.';«,Od.')^- (^2, s^^u, I; A) (35) 
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with the following boundary conditions: 



X^(0, s; u, I; A) + Xj^ki^, s; u, l;X) = ro(s) G Ak{u, I) 



(36) 



while the point (x, y) is the cross point of the respective trajectories belonging to different 
bundles, i.e. 



The vanishing superposition (35) if defined on the bundle By{u',l') suggests that 
^'^!^,(d, s; n', A) should be related somehow to the SWF '^'^{d,h^^{s;u,l);u,l);u,l; X) de- 
fined on the bundle Bk{u,l) which the previous one is a reflection. In the next section this 
relation is established as a condition matching both the solutions. 

4.5 Matching solutions defined on two different bundles of a skeleton 

Let now ^'Jf^''^ {x,y;u,l; Xys be defined in the above way on the corresponding bun- 
dles of B and B"^. Suppose that the corresponding SWF's ^'^((i, s; u, /; A), a = ±, and 
^l,''^{d,s;u,l;X), a = ±, which construct them can be continued along the respective bun- 
dles on which they are defined. Note that the domains Dk{u,l) and D^{u,l) where these 
SWF's are initially defined are then extended to DBk{u,l) and DB^(u,l) respectively. We 
are then faced with the problem of matching the continued solutions with the others defined 
on the same skeletons B and B"^. 

Anticipating the results of App.B such a matching should be done as follows. 

Let Bk{u,l) be a reflection of the bundles Bk^{ui,h),Bk2{u2,l2), ■ ■ ■ ,Bk„{un,ln) satis- 
fying (24). Let ^l_^{d,s;ui,li;X),'^l^{d,s;u2,l2;X),...,'^l^{d,s;Un,ln;X) be defined and 
continued on the respective bundles Bk^{ui,li), Bk2{u2,h), ■ ■ ■ , Bk^(un,ln)- Let further 
^'^°^((i, s; u, I; A) be defined on the bundle B^{u, I) while ^%{d, s; u, I; A) on the bundle Bk{u, I) 
being both related by the boundary condition (36). 

Then we make the following identification: 



r = [x,y] = rfe(di,si;n,Z) = ro(si) -|- di(si; n, Z) = 

rA;fc('^2, 82] U, I) = ro(s2) + d2(s2; U, I) 

rk{d, s; u, I) G Bk{u, I), rA;k{d, s; u, I) G B^{u, I) 



(37) 



^A%id,h{s;uj,ljy,u,l;X) = %.{D{s;uj,lj) - d, s;uj,lj; X) 
ro{h{s; Uj , Ij)) G Ak{u,l) Ci BL{uj,lj) 
ro(s) G L{uj,lj) 

ro{h{s;uj,lj)) = ro{s) + 'D{s;uj,lj), j = l,...n 



(38) 



Rewritten in terms of the x-coefficients eq.(38) gives: 



XA^ki^^^ Hs; Uj,lj)]u, I; A) = 




2 

X%{D{s\Uj,lj) - d,s-Uj,lj;X) 




ro{h{s;uj,lj)) G Ak(u,l) r\BL{uj,lj) 
ro(s) G L{uj,lj) 

ro{h{s;uj,lj)) = ro(s) + 'D{s;uj,lj), j = l,...n 



(39) 
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As previously 5{uj,lj) in the above formula is s-independent (see App.C). Due to that 
and due to the properties (20) and (21) the rhs of (39) satisfies (11) as it should. 
Putting d = in (39) and taking into account (36) we get: 

Xk{^,h{s;uj,lj);u,l;\) = -x]^ki^,h{s\Uj ,lj)\u,l;\) = 
dh{s; Uj, Ij 



ds 



Xk,iD{s;uj,lj),s;uj,lj-X) = 

_^^gaUp5K.,;,)^^cont(0^ /i(,; n„ /,); /; A) 

ro(/i(s; Uj,lj)) G A^iu, I) n BL{uj, Ij) 
ro(s) G L{uj,lj) 

ro{h(s;uj,lj)) = ro{s) + U{s;uj,lj), j = l,...n (40) 
Comparing (40) with (33) it is easily seen that ^"^{d, s;u,l; X) has to coincide with 

4.6 SWF's defined on a bundle skeleton 

Let r = (x,y) be a fixed point of the billiards. Let D{x,y) denote a set of all D{u,l), 
B{u,l) G B which contain this point and DA{x,y) the respective set of D'^{u,l), B^{u,l) G 
B"^. Let further DB{x,y) denote a set of all DB{u,l), B{u,l) G B which contain the point 
{x,y) and DB^{x,y) the respective set of DB^{u, I), B^{u,l) G B^. 

SWF's ^'g''^(x, y, A) vanishing on the billiards boundary can now be defined on B and 
B"^ as follows: 

*B^"(x,y,A) = 

^"{diu, I), siu, 0; u, l;X)+ Yl ^TKu, I), s(n, /); u, I; A) = 

D{u,l)eD{x,y) D^{u,l)eDA{x,y) 

Y *"(d(n,Z),s(u,Z);u,Z;A) = ^ ^J''{d{u,l),s{u,iy,u,l;X) 

DB{u,l)eDB(x,y) DB^{u,l)eDB^{x,y) 

r = ro{s{u,l))+d{u,l) (41) 

The SWF's ^''{d, s;u,l; X) defined on bundles of B and ^^(d, s; n, /; A) defined on re- 
spective bundles of B"^ are related with each other by the boundary conditions (36) and by 
matching conditions (38)- (40). 

'/^'^{x,y,X) aii^ -r^ 

It is clear that the conditions (40) have to determine also the ^-factors x^(s;n, /;A) 
for all the bundles Bk{u,l) which are the "initial" conditions for both Xfe (f^) -s; ^) and 
XA;kid,s;u,l;X), i.e. Xkis;u,l;X) = Xki^, s;u,l; X) = -XA;ki^,s;u,l; X). Nevertheless these 
conditions cannot be given arbitrarily. Just opposite all Xk{s;u,l; X) have to satisfy (40) in 
selfconsistcnt way. It is not easy to predict how this can be realized. 

Certainly one has to take into account that x^(s; u, I; A) being given at the point ro(s(n, /)) 
of the bundle Bi^(u,l) are distributed along the boundary dB by subsequent reflections of 
the initial ray r = r(j{s{u,l)) + d{u,l) which carries on Xki^: s;u,l; X) according the for- 
mulae (31), (33) and (40). However the ray come back eventually to the segment ^^(u,/) 
defining Xfc(s;'u,^;A) in some new point ro(s'(n, /)) of this segment. Naturally a value of 
x'^''^^{D{u,l), s' {u,l);u' ,1'; X) at this point has to coincide with Xfc(s'(^*! 0; ^! defined 



The solutions ^^'~^{x, y, A) and ' (x, y, A) coincide if and only if B = B"^. 
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already at this point. This is just one of the selfconsistency conditions contained in (40) (see 
"the last quantization condition" (49) below). 

It seems difficult however to find a general systematic way of exploring the system (40) 
to get solutions satisfying it. Nevertheless one can always think at least minimally trying to 
guess reasonable solutions to this system. In the next sections we will try to use both the 
approaches for some particular cases of billiards. 

The formulae (39) and (40) defines the conditions which the SWF's ^^(ci, h{s; Uj,lj);u, I; A) 
should satisfy when bouncing from the billiards boundary. Nevertheless this condition should 
be specified additionally with respect to its factors. Namely, taking their large A-limit we 
get: 



Xk,oih{s;uj,lj);u,l) 



dh{s; Uj,lj, 



1 

" 2 

aiXpS{uj a 



Xk,r+l{Ks\Uj,lj)\U,l) = -S 



ds 

dh{s; Uj, Ij) 



ds 



1 



/ (jj^ r-U(s;Uj,ij) / 1 / _i \ 

2 Er-i+iXkj,l{o., s; Uj,lj) da 
1=0 ) } 

r = 0,l,2,..., (42) 

The above equations should be satisfied on each bundle Bk{u, I) of the skeleton B. 

The first of the equations (42) should determine the classical quantities, namely the 
skeleton B and the " classical" energy Eq = ^p^ and by them define the JWKB approximation 
of the SWF's. Namely: 

*i^^^^^^(x,y,A) = Y ^'™^''id{u,l),siu,l)-u,l) = 

DB{u,l)£D(x,y) 

J2 J'Hd{u,l),s{u,l))e n ^ J" )xUd{n,l),s{u,l);u,l) (43) 

DB{u,l)eD{x,y) 

The remaining equations determine quantum corrections to the "classical" ones involved 
in (43). 

However it is easy to note that for the selfconsistency of the equations (42) it is necessary 
for the exponent e^'^^P^i'^j 'h) to be independent of A, i.e. we have to have on each bundle 
Bk{u,l) of B: 

Xp6k{u,l) = (pkiuj) 

Bk{u, I) CB (44) 

where Sk{u,l) is given by (33) and (f)k{u,l) is a A-independent constant. 

The equations (44) have to define both the skeleton B and the energy Eq = i^p^- 
Taking into account the last conclusions we get the following final set of the recurrent 

quantization conditions: 

Xp6k{u,l) = (pk{u,l) 
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Xl,r+iihis;uj,ljy,u,l) = _<5e'^^<^'=. 



dh{s; Uj,lj 



ds 



dh{s; Uj, Ij) 




ds 





pD{s;uj,lj) / I / I \ ^ \ \ 

r = 0,l,2,.., (45) 

together with: 

Xlfi{d,s;u,l) = Xk,o{s;u,l) 
Xk,r+i{d, s; u, I) = Xk,r+i{s; u, I) + 

^lo if^'^ {-^'^XkA^^ ^' O) + 2 X] Er-i+iXk,m{a, s; u, Z) j da 

r = 0,l,2,..., (46) 



4.7 Continuity of SWF's built on a skeleton 

A skeleton on which SWF ^"-^{x, y, A) is built in the way described above does not ensure 
automatically that this SWF will be defined continuously on it. It is certainly continuous 
inside each bundle of the skeleton but it can appear to be discontinuous if a point (x, y) 
crosses bundle boundaries. In general the following circumstances can accompany such a 
crossing: 

1. the crossed boundary of DB^iu^l) coincides with a boundary of DB^iiu + l^V) of a 
neighboring bundle Bj^i{u + Z, V); 

2. there is no such respective neighboring bundle. 

In the first of the above cases the corresponding ^'^((i, s; u, Z; A), {x(d, s),y{d, s)) G 
DBk{u,l), and ^'l:,{d,s;u + l,l';X), {x{d, s),y{d, s)) G DB^'in + 1,1'), defined by (35) have 
be identified on the common boundary, i.e. 

*^(d, s; u, I; A) = ^'^,((1, s;u + l, l'; A) 
{x{d,s),y{d,s)) edDBk{u,l)ndDBk'{u + l,l') (47) 

In the second case however the corresponding SWF's "^^{d, s;u,l; X), {x{d, s),y{d, s)) G 
DBk{u, I) have to vanish on such a boundary of DBk{u, I), i.e. 

^^{d,s;u,l-\) = 
{x{d,s),y{d,s))£dDBk{u,l) (48) 

The last condition though necessary seems to be a little bit arbitrary. However we should 
remember that our calculations are performed in the semiclassical regime, i.e. in the classi- 
cally allowed regions (bundles) outside which the semiclassical wave functions cannot exist. 
Physically this means obviously that outside each bundle a corresponding piece of the wave 
function represented on the bundle by its respective semiclassical approximation has to vanish 
exponentially when moving away from the bundle. 
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It will appear in the next sections that both the possibilities (47) and (48) will be met 
and will have to be applied to solve the problem of the semiclassical quantization. 

By the definition (41) and by the recurrent equations (45)-(46) and by the conditions 
(47)- (48) we have formulated the quantization procedure which should determine skeletons, 
the SWF's ^'g '"^(x, y, A) and the corresponding energy levels. In the next two sections we 
shall apply this procedure to the simplest well known cases of billiards, i.e. to the circular 
and the rectangular ones. 

Let us note finally that if B 7^ B^ then energy levels corresponding to the skeleton B have 
to be degenerate. This conclusion follows easily from the form of the quantization conditions 
(45)-(46) and (47)-(48) showing that the complex conjugations of \I'g '"^(x, y. A) satisfy also 
these conditions with the same semiclassical energy E. The two corresponding solutions are 
of course ^^'^{x,y,\). 

4.8 Finite and infinite bundle structures of skeletons. The last quantiza- 
tion condition 

For a given billiards there can be skeletons with a finite number of bundles as well as with 
an infinite one. The semiclassical quantization procedure described in the previous sections 
seems to be easily applied to the finite bundle number skeletons. Namely in such a case 
following a trajectory starting from a bundle Bk{u,l) wc have to approach the same bundle 
after a finite number of bounces. The corresponding semiclassical wave fTinction propagated 
by the skeleton has therefore to come back to its initial form achieving again the initial 
bundle. This condition closes the process of quantization formulated in the previous sections. 
The respective conditions are of course the following: 

(-<5)-exp[ <Pk{u,l)\=l 

xT"\D,s-u,h\) = xl{s;u,l-\) (49) 

where n is a number of bounces and D is the total distance passed by the billiards ball along 
the investigated trajectory. 

The cases of skeletons with an infinite number of bundles are much more difficult for 
investigations. Such skeletons should be typical for chaotic billiards. According to its defini- 
tion a bundles Bk{u, I) can bifurcate after the reflection by the billiards boundary into many 
different subbundles, i.e. parts of other bundles having their beginnings also partly on the 
arc BAk{u, I). In fact a general behaviour of a skeleton in such chaotic cases should not differ 
essentially by its chaotic complexity from a chaotic trajectory reminding however rather a 
gigantic road-knot with inflnitely many viaducts spanning the billiards boundary on which 
the billiards ball moves. It is obvious that if they exist their identification seems to be not 
an easy task. 

Nevertheless the rule (49) can appear to be useful also even in such cases. This is because 
a ray beginning with a bundle Bk{u, I) can come back to it even arbitrarily close to its initial 
starting point on Ak{u,l) (according to the Poincare theorem) not repeating its way. But 
this is enough for writing the "last quantization condition" (49) where the sum goes now over 
all bundles of the skeleton passed by the ray. 

In the next sections we shall focus on the finite number cases of bundles in skeletons not 
avoiding however billiards with chaotic motions such as the Bunimovich one. We shall discuss 
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Figure 2: The two bundle S(0,27r;a) and B^(0,27r; vr — a), < a, of the skeleton B in the 
circular billiards 



also in sec. 8 a possibility of explaining by the skeleton idea phenomena of scars. These are 
just phenomena which cannot certainly be described by skeletons with finite structures. 

5 The circular billiards 

In this section we would like to apply the procedure of constructing the SWF's on the 
classical trajectories according to the last section to the circular billiards Fig. 2 and to quantize 
with these rules the energy of the circular billiard ball. 

First we have to construct a skeleton B. It follows from the considerations of App.B that 
any of it contains only a single bundle -6(0, 27r; a) defined on the full circle L(0, 27r) rays of 
which make all the same angle a, < q < ^vr, with the tangents to the circle. Its partner 
B"^ contains also a single bundle -8(0, 27r; vr — a) the rays of which make the angle vr — a with 
the tangents to the circle. Therefore in the case a = both the bundles coincide. 

Motions along such rays conserve of course the angular momentum pcosa so that both 
the bundles respect the cylindrical symmetry of the billiards. 

For these two bundles the corresponding SWF's are the following (see Fig. 2): 

^^{d,s,X) = ((i-sina)-5e±*^P('='±^™^")x^(rf,s,A) 

^i{d,s,X) = {d-sma)-h^'^P'^'^^'^°'''^xi{d,s,X) (50) 

The above SWF's are defined in the ring < d < sin a, < s < 27r of the circular 
billiards. 
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Since both the bundles are defined in the ring mentioned the SWF's have to satisfy the 
fohowing uniqueness conditions: 

^±(d,s,A) = *=t(d,s + 27r,A) 

*±(d,s,A) = *±(d,s + 27r,A) (51) 

which lead to 

Apcosa = — Apcos(7r — a) = m, m = 0, ±1, ±2, ... (52) 
and means the angular momentum quantization together with 

X^{cLs + 2TT,X)=x^{d,s,X) 

xiid,s + 2Tr,X) = xi{d,s,\) (53) 

We can demand from the SWF's ^^{d, s, A) and ^^{d, s, A) to be eigenfunctions ^±m{d, s, A) 
(= "^^{d, s, X)) and A:±m{d, s, X) (= "^^{d,s,X)) of the respective angular momenta ±mh, 
m = 0,1,2,..., where m = Xpcosa. Then it follows easily from the forms of the SWF's 
considered that for the corresponding x±m{d, s, A) and XA;±m{d, s, A) we have to have: 

dx±m{d,s,X) ^ dxA;±m{d, s, A) ^ ^ 
ds ds 

i.e. both x±m{d, s, A) and XA;±m{d, s, A) are independent of s and therefore they are constant 
on the circle boundary, i.e. for d = 0. 

Putting therefore Xto(0, s. A) = 1 and assuming further the following boundary conditions: 

X±m(0,S,A) = -XA;±m(0, S, A) = 1 (55) 

we get for the appropriate ^'5f^(a;, y. A) vanishing on the circle boundary and being the 
angular momentum eigenfunction: 

*Sn(x, y, A) = '^±mid, Si, A) + *A;±m(c?, S2, A) 

m = 0,l,2,... 

(x, y) = (cos si, sin si) + d{— sin(a + si), cos(a + si)) = 

(cos S2, sin S2) — d{— sin(— a + S2), cos(— a + S2)) (56) 

Making further the identification 

*A;±m(d, s, A) = 'I'±^*(2 sin a- 2a, A), sina>d>0 (57) 

where ^'j?^*((i, s. A) denote '^±m{d, s, X) continued by the caustic r = cos a we get for 
*5f^(a;,y,A) = *5?^(r cos$,r sin$, A) (see Fig.2): 

*5f^(rcos$,rsin$,A) = ^'±^(ci, s. A) + ^'A;±^(d, 2$ - s. A) = 

^±m{d, s, A) + ^^\2 sin a - 2$ - s - 2a, A) = 

(sina-d)-^ (-ie^'^P''^'"''x±mid,s,X) ± 



gTiAp(d-2sina)±jm(2*-s-2a)^cont^2sina - d, 2$ - S - 2q, A)j = 

m = 0,1,2,.. 



d = sin a — \/ — cos^ a 



sin($ — s) = (sin a — \J r"^ — cos^ a) 

sina>d>0 (58) 
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The identification conditions (57) give: 

XA-±m{d, S, A) = e±2i(Apsina-ma-|)^ca.t(2 sma-d,s- 2a, A) 

sinQ;>d>0 (59) 

Applying now (49) and (55) we get: 

p±2j(Apsina— ma+^) 

xS^*(2sina,s -2a,A) = x±m(0,s,A) = 1 (60) 

so that 

/■2sina / , , \ 

-l+lX±m;l{0', s) \ da = 
■'^ \ l=Q J 

A; = 0,1,... (61) 

with 

X±m-fi{d,s) = 1 

X±m;k+l{d, S) = (^J^^[j~^X±m;k{(^,s)'j + 2^^ Ek-l+lX±md{(^, s)^ da 

k = 0,1,2,... (62) 

The integrations in the formulae (61)-(62) go above the singular point d = sin a in the 

d-plane for the plus sign and below this point for the minus one. 

A simple conclusion from the quantization conditions (61)-(62) is that: 

{X+m;k{d, s))* = X-m;k{d, s) 

k = 0,1,2,... (63) 

so that 

{n^^{x,y,X)y = ^^_Ux,y,\) (64) 

i.e. the semiclassical energies E in the circular billiards are degenerate with respect to the 
sign of the angular momentum. 

Of course the above conclusions are well known as exact for the circular billiards. 

Therefore in the equations (61)- (62) we can choose only the plus sign in the respective 
calculations to get: 



/ m TT 

XfEo.,mr - ^ - marccos = ttt - - 

m 

COS (^fYiv 



kmr 



m = 0,l,2,..., r = l,2,... (65) 



for the "classical" energy Eq and the angle a. 
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It follows from (65) that 







COS CCor- 



m 



r- i 



cosa^r = F( -T^) 

m 

m = l,2,..., r = l,2,... (66) 

where F~^{x) = {^/l — x'^ — x aiccos x)x~^ , < x < 1. 

The next to zero order terms of the corresponding semiclassical expansions can be obtained 
using the results of App.D from which it follows (the formulae (77)- (79)) that the operator 
D{t, s) = Ja • A • J~ 2 in (68) is bilinear in ^ and ^ with coefficients completely independent 
of s and it acts in (68) also on the s-independent quantities. In such cases its action is reduced 
to (for t ^ X = pt): 



red ' dx' 



(cos^ a ^\ 9^ 2 cos^ ad 5 cos^ a 11 

where J = sin a — x. 

Taking this into account in the second of the eqs. (61) and putting there A; = we get: 



2 (•2 sin amr , ^ 1/5 \ 

El-mr = --r-. / I J2 • A • J"2 ) ■ldx= — - COt^ amr - 1 

4smQ;^r ^0 ^ ^^ed 8 sm amr \6 / 



(68) 



where the integration in the x-plane went over upper half-circle with the center at a; = sin a^r 
in the clockwise direction. 

Similarly for Xmr;iid) we get: 

i (5 cos^ amr I I \ i 




J^r- 12 sin^ a 




2p JO y4 Jmr ^ "^rnr j P 

5 cos^ amr 5 cos^ a. 



(69) 



where Jmr = sin amr — d. 

The higher order terms of the semiclassical expansions for EmrW and Xmr{d,X) can be 
obtain analogously using the recurrent equations (61)-(62). 

It is to be noticed that the above calculations form a new algorithm for the semiclassical 
approximation method. 



6 The rectangular billiards 

Consider now the rectangular billiards shown in Fig. 3. This billiards is the canonical 
example of the energy quantization problem because of its easiness to be solved by the variable 
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Figure 3: The rectangular billiards with eight ray bundles 



separation method. According to Fig. 3 the well known solution to the problem is given by 
the following two equations: 

sm{Xpxa) = 

sm{Xpyb) = (70) 
with Px+Py=p'^ = as the result of the following form of the energy eigenfunction: 

^{x,y) = C sm{Xpxx) sm{Xpyy) = 

__ (^PxX+\pyy _|_ ^-Xp^x-Xpyy _ ^Xp^x-Xpyy _ ^-Xp^x+Xpyy^ ^j-^^ 

Of course one can always put px = p cos a, py = psina where a, < a < ^vr, is the angle 
the momentum p is inclined to the x-axis. 

Let us note that the cases a = 0, ^vr are excluded by the solutions (71). 



6.1 The cases < a < 

Let us consider this rectangular billiard example by the method described in the previous 
sections. Since the absolute values of the momentum components Px,Py are the integrals of 
the classical motion inside the billiards respecting elastic law of bouncing then all bundles 
which should be taken into account are defined by a single angle a, < a < ^vr, which are 
made by the rays of the bundle Bi = B{0, a; a) with the x-axis, see Fig.3. 

Choosing the case of the angle a shown in Fig.3 the remaining seven bundles of the 
skeleton B shown in Fig.3 are B2 = -6(0, a; vr — a), B3 = B{a,b;7r — a), B^ = B{a,b;TT + 
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a), = B{a + b,a;iT + a), Bq = B{a + b, a; 2tt — a), B-j = B{2a + b, b; 2tt — a), Bs = 
B{2a + b,b;a), i.e. the parameter s introduced in sec. 3 is counted anticlockwise starting 
from the point (0, 0) of Fig.3 (and having negative value if measured clockwise). The bundles 
B2k-i, B2k, k = 1, ...,4, are defined on the respective sides Lk, k = 1, ...,4, of the billiards, 
i.e. on Li = L(0, a), L2 = L{a, b), L3 = L{a + b, a), L4 = L(2a + b, b). 

The skeleton coincides exactly with B in the case of the rectangular billiards. 

Let us note that a number of bundles in the skeletons is obviously independent of a choice 
of a, i.e. it is always equal to eight. 

Note also that unlike the circular billiard bundles the rectangular bundles do not have 
their caustics inside the billiards, but rather in the infinities. 

The semiclassical wave function ^'"^(x, y) constructed according the rules discussed earlier 
should be the sum at most of the eight SWF's 'i'f{d, s, X), i = 1, ...,8, which arc constructed 
on the respective bundles Bi, i = 1, 8, of the skeleton B. But as it follows from Fig. 4 only 
four of them interfere in each point of the rectangular billiards, i.e. we have: 



where Q:(.s,j) arc the angles the rays in Fig. 4 make with the x-axis and (xo(si), yo(si)) are the 
points of the rectangular boundary from which the rays start. The parameters Sj as previously 
measure distances of the boundary points {xo(si) , yo(si)) from the point (0, 0) in anticlockwise 
direction along the boundary. The corresponding Jamaicans are J{t, Sj) = -\/— sina(si), i.e. 
are constant but discontinues. Therefore they will be included into the x-factors contained 
in the SWF's. 

Let us note that y) can be discontinue each time the point (x,y) crosses the thin 

black lines in Fig.3. being the boundary of different bundles. Let us enumerate these lines 
by li, i = 1,...,4, starting from the line emerging from the point (0,0) and continuing 

anticlockwise. 

It is then easy to see that to ensure the continuity of \l'"*'(x,?/) in each point in the 
billiards it is enough to identify the respective x-coefficients on the respective lines defining 
them elsewhere in the billiards as continuous functions of s and d. We get such a continuity 
putting: 



8 4 



^'-(x, y) = Y^ *+(d,, s^, A) = ^ (d„ s„ A) 



i=l i=l 



X = Xo{si) + di COS a{si) 
y = yo{si) + disma{si) 
i = l,...,4 



(72) 



xfid,0,X)=xt{d,0,^) 
X3{d,a,X) = xt{d,a,X) 
xf{d, a + b,X) = xi{d, a + b,X) 
X^{d, 2a + b, X) = xi{d, 2a + b,X) 



(73) 



Of course *"*(x,y) vanishes on the rectangular boundary, i.e. we have to have: 



^«^(0,j/) = *''"(o,y) = ^""(x,0) = *"*(x,6) = 

0<x<a, <y <b 



(74) 
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Figure 4: Four solutions "^f{di,Si), i = 1,...,4, meeting at the point 



The conditions (74) expressed in terms of the SWF's ^f{d,s,\), i = 1,...,8, take the 
following detailed forms: 



*+_,(0,s,A) + ^2-,(0,s,A) =0 
*+(0,s,A) + ^'2,_i(0'^'A) =0 
s G Lk, k = l,...,4 



(75) 



When identifying ^'^^(0,5, A), i = 1,...,8, with the respective "^^{dj{s), fj{s), X), j = 
1,...,8, it should be noticed that for each considered trajectory family there is no caustic 
inside the rectangular billiards so that each integral Jj^pxdx + Pydy along the close loop 
K lying inside the billiards or on its boundary has to vanish. As a consequence of this 
the identifications we are talking about are then reduced just to the identifications of the 
corresponding x-factors, i.e the possible phase factors are factored out. Therefore we get 
successively for the conditions (75) to be satisfied: 

1. s > 0, ro(s) G /, -stauQ G XII, 2a + b- 6cot a - ro(s) € VII 



Xi (0,s, A) = xt 



Vsina 



2a + b — b cot a — s,X 



and hence and from (75) 



X2 (0,s, A) = Xt 



Xi(0,s, A) +X7 



/ s 



, cos a 



-s tan a, A 



(76) 



cos a 



-s tan a, A = 
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X2(0,s,X) + X4(^^,'^a + b-bcota-s,x] =0 (77) 
Vsma / 

In the last two equations the plus signs at x^, j = 1, 8, have been dropped and this 
convention will be kept in further equations if not leading to misunderstandings. 

2. s > 0, ro(s) G //, -Stan a G X + XI, a+ (a - s)tana G V + VI) 

Xi{0,s,\) = xt {- — -,a + {a - s)tana,x) 
\ cos a ) 

X2 (0,s,A) = X7 (^^>-«tana,A) (78) 

and hence and from (75) 

Xi(0,s,A) + X7 (^—, -Stan a, A) = 
V cos a J 

X2(0,s,A) + X4 [- — ^,a + (a- s)tanQ;,A) = (79) 

\ COS Oi J 

3. ro(s) G ///, s > 0, (-6- s + 6cota G IX, a + (a - s)tana G IV 

Xr(0,s,A) = X4 (- — ^,a + (a-s)tana,A) 

\ COS J 

X2 (0, s, A) = xt -b-s + bcot a, x) (80) 

V sm a / 



and hence 



Xi(0, s, A) + X6 f^— j-^*- s + 6cota, A^ =0 
V sm a J 

X2(0, s, A) + X4 ( - — -,a + {a - s) tana, A ) = (81) 
V cos a ) 



4. s > 0, ro(s) G /y + y + y/, a - (s - a) cot a G /I + ///, 
a + 6 + (a + 6 - s) cot a G Vll + F/// 



Xs (0,s,A) = Xe f^-T^ — ^,a + 6 + (a + 6- s)cotQ;,A) 
V sma / 



X4 (0, s, A) = xM a - (s - a) cot a, A ) (82) 

' sma ' 



and hence 



X3(0, s)e-^^?"^'^°^° + XI f - — - is - a) cot a, A^ e'^^"*^"^" = 

V sm a / 

X4(0, s, \y-^^vaco^<^ " + + 6 + (a + 6 - s) cot a, A^ e^^^"*""^" = (83) 

V sm a / 
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5. s > 0, ro(s) G VII, 2a + b — 6cot a — ro(s) el, a + b— {s — a — b) tana G 

b 



X5(0,s,A) = x^ ( 



sin a 



-, 2a + 6 — 6 cot a — s,X 



Xe (0,s,A) =X3 (^^^^^^,a + 6- (s-a-6)tana,A^ (84) 



and hence 



X6(0, s, \)e-^^P^'^'' + xi( -^,2a + b - bcot a - s, \) e^^f^^^"" = 

\sina / 

X5(0, s, A)e-^^f^«'°" + X3 f ^^^^—^,a + b - (s - a - b)tana, x) e'^^^'"''"' = (85) 

V COS a / 



6. s < 0, ro(s) G VIII, -6- (6 + s) tan a e X + XI + XII, 
a + b - (a + b + s)t&-n.a e IV + V 



X3 (0,s) = x^ (— — -,-b - {b + s)ta-Q.a,\\ 

\ COS Oi J 

X4 (0, s) = xt — ~~ ' c^ + b— {a + b + s) tan a, (86) 



and hence 



X4(0,s,A)e-^^^'^^^°" + Xi (6 + s)tana,A) g+Upbsina ^ q 

V cos a / 

X3(0,s,A)e-^^f^^^°" + X2 ("^^^^^,0 + 6- (a + 6 + s)tana,A^ gUpfesina ^ q /g7^ 

V cos a / 

7. s < 0, ro(s) G /X, -b- (6 + s)tanQ; G X, 6cot a - 6 - ro(s) G //J 

X5 (O'S) = Xs {— — (6 + s) tan a, a) 
\ cos a / 

X^{^,s) = xt{^,hcoia-b-s,\\ (88) 
V sm a J 

and hence 

X6(0, s, X)e-'^P''^^'' + X8 f Z^^, -b-{b + s) tana, A^ 6+^^^'^^^"" = 

V cos a / 

X5{0,s,\)e-'^P^'^'' + X2 ( -^,bcota-b-s,x] g^^f^'^i'^" = (89) 

Vsma / 

8. s <0, ro(s) G X + XI + XII, -scot a G / + //, 
- 6 - (s + 6) cot a G y/// + /X 

X8 (0,s,A) = xt (s + 6)cota,A') 

Xf(0,s,A) =x^ (^,-scota,A') (90) 
\ sin a / 
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-b-(s+b)cota 




Figure 5: The solution s) being carried by the successive bundles Bg 

^ B. ^B-r^B, 



Ba 



Bi 



and hence 



X7(0,s,A) +X2 

X8(0,s,A) +X3 f^T-^,-fe 
V sm a 



-, —s cot a, A 







sma 

— {s + b) cot a, A ) =0 



(91) 



Solving the last set of the quantization conditions we should get the solutions for the 
energy levels and for the corresponding ^""^(x,!/). 

However the above conditions do not provide us with the solutions for Xj(0, s,A), i = 
1,...,8, in some algebraic forms. Rather they show how a solution defined on some bundle 
propagates along its rays to be transformed by bouncing off the rectangular boundary into 
another bundle and finally after a finite number of such bouncings achieving its mother bundle 
and repeating this process infinitely. 

To get a knowledge what really happens by such bouncings consider for example the 
green color ray along which the solution ^'g'(t,s) propagates, shown in Fig. 5 emerging from 
the point s, ro(s) G X, of the rectangular boundary. This particular ray bounces into the 
blue-light one and next again into the green one. Taking into account first the conditions 
(85) and next (81) we get: 



Xi(0, (s + 26) cot a, A) = -xe 



sm a 



-6 — (s + b) cot a, A 



cont 



/ 2b + s 
V sin a 



(92) 
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and continuing the bounces to achieve the green ray again we get: 

Xi(0,(s + 46)cota-2a,A) =X8°"*(^(s),s,A)e^^^J'''^^'^"+2^^^'«^°^", ro(s) G X (93) 

where D(s) is the total distance the green ray has passed from the point s to the point 
(s + 46) cot a — 2a bouncing muhiply from the boundary. 

But taking into account the first of the identification (73) and the identification (49) we 
have to have: 

Xi(0, (s + 26) cot a, A) = Xg""* f s, x) 

XI (0, {s + 46) cot a - 2a, A) = xr*(^(s), A) 

ro(s) G X (94) 



so that 



and consequently: 



g+2iAp6sina j 

g4jAp6sina+2iApacosa j (95) 



Ap6 sin a = m7r 
Apa cos a = titt 
m,n = 1,2,... (96) 

Taking the limit A ^ oo in (94) we get: 

Xi,o((s + 26) cot a) = Xi,o((s + 46) cot a - 2a) = X8,o('S) (97) 

If tana = is not rational none of the arrival points can be repeated, i.e. the corre- 
sponding trajectory is not closed (periodic), and the (infinite) set of such points is dense on 
the side L\ of the billiards. Therefore the equations (97) show that Xi,o(5) has to be not only 
a constant of motion but also s-independent since this equation can be written in infinitely 
many points of L\ densely distributed on it. 

The last conclusion is valid also of course for the remaining Xfc,o(s), fc = 1, 8, and as it 
can be easily concluded from (14) and App.D all the coefficients Xk{d^ s, A), k = 1, 8, are 
then also constant. Putting therefore xiid:S,X) = 1 we get by (73) and (75): 

Xk{d,s,X) = l, A; = 1,4, 5, 8 
Xk{d,s,X) = -l, A; = 2, 3, 6, 7 (98) 

Choosing therefore the point {x,y) of Fig.4 for the SWF (72) we get: 

*«^(X, y) = *+(di, Si) + *+(d2, S2) + (4, S3) + *+(d4, S4) = 

4 

1=1 

where ti and fi{si), I = 1, ...,4 should be calculated from the relations (see Fig.5): 

X = xo{si) + pti cos a{si) 
y = yo{si)+pti sma{si) 

a{si) = a, -TT — a (100) 
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However making use of the independence of the phase integral JI^qq^ Pxdx + Pydy of the 
integration paths we get for the particular terms in the sum in (99): 

giAp2<i+jAp/i(si) _ ^iX{-pxX+Pyy) 

^iXp^t2+i\pf2{s2) _ ^i\{-pxx-pyy) 
^iXp'^ts+iXpfii-ss) _ ^i\(pxX-pyy) 

^iXp^ti+iXpfiisi) _ ^iXiPxX+Pyv) (101) 

and finally: 

*"*(x,y) = —e^^^~'P'^^'^'Pyy^ + e^K-'P'^^'PyV) _ ^iKP^^^-PyV) _|_ gi>^{pxx+pyy) _ 

-4sm{pxx)sm{pyy) (102) 

reproducing in this way the exact result (71). 

The cases when tana = ^| is rational are possible only if | is rational. For maintaining 

the results for the irrational | one can argue relying on a continuity of all the investigated 
quantities considered as functions of a, b, since each point (a, b) with rational value of | is 
densely surrounded by the ones for which | is irrational. 

6.2 The cases a — 0, ^tt - the bouncing mode skeletons 

It is surprising that the cases a = 0, ^tt are not allowed by the representation (71) of the 
SWF leading to the totally vanishing solutions while there are still the bundle skeletons with 
these angles on which non vanishing identically SWF's can be constructed. Therefore we 
should get some new knowledge about the semiclasscical method developed here considering 
these cases known as the bouncing modes. 

Consider the case a = ^ir. The second case will be analogous. 

There are two bundles in this case which the skeleton B is consisting of. The one Bi with 
its rays directed up and starting from the side Li and the second with rays directed down 
starting from the side L3 (Fig. 6). The skeleton B"^ is identical with B. 

For these particular cases of bundles rays for both the bundles will be positioned by the 
same parameter s measuring a distance of a ray from the y-axis along the corresponding sides 
Li and L3. Therefore for the corresponding SWF's we get: 

^i{b - d, s, A) = e^'^P^'-'^hfib - d, s, A) 

0<d<b, 0<s<a (103) 

For the SWF ■^"■%x,y, A) we have: 

*"^(x, y, A) = ^tid, s, A) + *+(6 - d, s, A) = *r(6 - d, s, A) + ^^{d, s, A) 

{x,y) = {s,d) (104) 

together with the following identifications: 

*^(6-d,s,A) = '^f{d,s,\) 
^'+(0,s,A) = -*r(0,s,A) 

^'+(0,s,A) = -*3(0,s,A) (105) 
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d=y 



b-d 



s=x 
B, 



Figure 6: The two bouncing mode bundles Bi and of the vertical skeleton in the rectan- 
gular billiards 



so that: 

xt{0,s,\) = -Xl{0,s,\) 

xt{0,s,X) = -X3{0,s,X) (106) 
It follows from the last equations that: 

X^(0, s, A) = -X3 (0, s, X) = -e^^P'xtib, s, A) (107) 

so that: 

xt{b,s,X) = -e'^'''xti'2b,s,X) (108) 

and further: 

(0, s, A) = -xr(0, s, A) = -e'^^'xtib, s, A) = e^^^^^x^ (26, s, X) = e^^^^VlO, s, A) (109) 

since XiC^b, s, A) = x^™"*(0, s, A) = Xi{0, s, A). 
Therefore: 

and 

^«^(x, y, A) = e'^P^xtid, A) + e-'^P^e'^^'xt (b - d, s, A) = 

e^^P'^xtid, s, A) - e-^^'^x^ (26 - d, s, A) (111) 
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We have of course ^'"*(a:,0, A) = ^"^(a:,6, A) = by construction. But we have to have 
also ^"^(0, ?/,A) = \L'"*(a,y,A) = 0, i.e. we have to have: 



e 



iXpd 



xtid,0,X)-e~'^P''xti'^b-d,0,X) = 



e'^P\tid, a, A) - e-'^^'xti^b -d,a,\)=0 (112) 
Therefore from (112) in the hmit A — > oo we have: 

e^'^'xto(«)-e-*'"'xto(«) = (113) 

so that 

xto(O) = xto(«) = (114) 
Next let us invoke the second of the equations (14) to get in the considered case for d = 2b: 

xU2b, s) = xUO, s) + + 2^ixto(^) j (115) 



so that 



+ 2Eixt,o{s) = (116) 



since Xi,i{'^b,s) = xli{0,s). 

The obvious solution of the last equation satisfying the boundary conditions (114) is: 

X^;o(s) = ^osin(v^s) 



V2^a = m7r, m = l,2,... (117) 

Coming back to the second of the equations (14) we can conclude that xtiid-i *) again is 
independent of d. Passing next to the third of the equations (14) and repeating arguments 
similar to those which led us to (115)-(116) we get the following explicit dependence of xt,i{^) 
on s: 

E2A0S 



xtiis) = Ai sin(V2^s) + Bi cos{V2Eis) + " cos{V2Eis) (118) 

The boundary conditions Xi^i(O) = Xi^i('^) = enforce however Bi = E2 = 0. 
Using again (14) and the inductive arguments we come to the conclusion that xt{d, s, A) 
is d-independent and coefficients of its semiclassical series have the form: 



Xlkis) = AksiniV^Eis), A; = 0,1,... (119) 
so is the form of Xi^(^) X) itself, i.e. 

X^-(s,A) = A(A)sin(v/2^s) (120) 
Therefore coming back to (111) we get: 

^"■'{x, y, A) = {e'^P'^ - e-'^P'^)xt{s, A) = 2i^(A) sin(Apd) sm{^/2E^s) (121) 



30 
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O 



Figure 7: An "arbitrary" broken rectangular billiards 



which again is the result of the previous way of the rectangular billiard energy quantization. 
But now the energy E is given by the (finite) semiclassical series: 



Therefore we get a surprising result that SWF's in the rectangular billiards can be built 
equivalently by the following two ways: 

1. on the skeletons which rays are inclined to the billiard sides with corresponding angles 
defined by the quantization conditions (95) - in this case the x-coefhcients of the SWF's 
are simply constant on the rectangular boundary; 

2. on one of the two skeletons which rays are perpendicular to one of the rectangular sides 
- in this case the corresponding x-coefficients vary along the sides perpendicular to the 
rays and vanish on the sides parallel to them. 

But a more important conclusion which follows from the results got in this section is 
that a relation between the form of skeletons and the SWF's and the corresponding energy 
spectrum which seemed to be suggested by the first part of the section is rather illusory since 
a full description of these quantities can be also obtained using a single skeleton only at least 
in the case of the rectangular billiards. 

6.3 Broken rectangular billiards 

By a broken rectangular billiards we understand the one on Fig. 7, i.e. with some number 
of rectangular bays and peninsulas. 




m, n = 1, 2, ... 



(122) 
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Figure 8: A single bay rectangular billiards 



To illustrate the way of energy quantization in such billiards we shall consider first the 
one with a single peninsula shown in Fig. 8. The corresponding procedures are analogous to 
the previous ones. One needs to construct additional four bundles relating with the peninsula 
sides. 

Although both the previous methods of the sections 5.1 and 5.2 are equivalent leading us 
to the same results the method of sec. 5. 2 seems however to be simpler and more instructive 
in applications to more complicated cases of billiards. 

We assume from the very beginning that all the sides of the billiards from Fig. 10, i.e. 
a = l,b,a',b' are commensurate. This assumption can always be satisfied even if the sides 
are expressed by irrational number just by respective approximations of the latter by rational 
ones with arbitrary accuracies. We proceed as follows. 



We construct two vertical skeletons. One for the rectangle OAEF and the second for 
ABCD. 

For the first skeleton according to (121) we get the following SWF: 



"^l^x, y, A) = Ai sm{Xpy) sm{^/2E^x) 
0<x <a', 0<y <b 
Xpb = mvr, m = 1,2, ... 



\/2E\(x = nvr, n = 1,2, ... 



(123) 



while for the second: 




(124) 
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To get (2^1 2/1 ^) foi" the total rectangular we have to match both the previous ones on 
the segment AD of Fig. 8. Since this matching has to be valid for < y' <b' then as it follows 
easily we have to have p = p' so that Ei = E[ because of (122), i.e. the vertical wave lengths 
of both the matched solutions have to be the same. Therefore we get: 



Ai sin( V2^a') = A2 sin V2^(l - a' 



Ai cos{y/2Ela') = -A2 cos (v^2^(l - a')) (125) 

so that ^/2E^ = Ztt, Z = 1, 2, and A2 = (-1)'+Mi. 

Therefore the procedure leads us to the following quantization conditions for the energy 



1 b Ay 

— = n — , — = m— ^ 

no 2 mo 2 

_ 27r A _ 



1 o 2E 27rV 1 1 \ tt"^ I o o m^ml 
2 + X2 



P _-„2, ^ , ^ _ " f ^2^2 I '"- '"-0 



m,n= 1,2,... (126) 

where A^., Aj, are the wave lengths of rays in the horizontal and vertical skeletons respectively 
and no, mo are the smallest integers satisfying = noa' and kob = mob' where lo, ko are also 
integers. 

The respective SWF's are the following: 

^as A) = I ^ '''' S 5 = ^ sin(nno7rx) sin (=y) {x, y) G D^r ^^^7) 

where denotes the domain of the x,y-planc occupied by the broken rectangular of Fig. 8. 

One can easily realize that the last results can be easily generalized to any broken rectan- 
gular billiards. A little bit surprising is that the formulae (126) for the energy and (127) for 
the wave functions remain unchanged for any such a billiards while a number of conditions 
the wave lengths A^; and A^ have to satisfy filling the vertical and horizontal skeletons by 
integer numbers of their halves is increasing respectively to numbers of bays and peninsulas 
forming the sides od such billiards. 



7 The simplest SWF's solutions for polygons and Bunimovich 
billiards 

The semiclassical quantization of the rectangular billiards in sec. 5. 2 shows that it is possi- 
ble to quantize semiclassically in the similar way some specific configurations of the SWF's in 
more complicated billiards also classically chaotic such as polygons or Bunimovich billiards. 
Such easy opportunities appear if billiards to be considered possess boundaries which allow us 
for easy constructions of corresponding skeletons. Some simple examples of such billiards are 
provided by a parallelogram, a trapezium, a pentagon shown in Fig. 9 or by the Bunimovich 
stadium - Fig.lO, and its simple generalizations - Fig.ll. 
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Figure 9: A parallelogram, a trapezium and a pentagon billiards with the bouncing mode 
skeletons 
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Figure 10: The Bunimovich billiards with the bouncing mode skeleton and with the scar 
generating bundles associated with the horizontal periodic orbit 
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The corresponding bundle skeletons are shown on the respective figures. 
The constructions of the SWF's corresponding to the skeletons shown on the figures 
coincide with those for the rectangular billiards and therefore these SWF's are the following: 



for the skeletons perpendicular to the x-axis and defined in the respective rectangles a x b 
of the figures while outside the skeletons y, A) = 0, i.e. those regions are classically 

"forbidden" and the corresponding wave functions vanish there exponentially if A ^ oo. . 

For the Bunimovich-like billiards of Fig. 11 it should be clear that bouncing modes can be 
excited independently in each of its rectangular part but also simultanuously in both parts 
tunning the modes respectively. 

Formulae like (128) have been suggested by Burq and Zworski [20] by completely different 
but mathematically rigorous arguments. 

8 Scars and periodic orbits 

It was shown in the previous section that the idea of the skeletons seems to be effective 
in solving some simple situations of quantum phenomena related semiclassically with the 
chaotic dynamics. Nevertheless obvious difficulties in effective constructions of skeletons in 
the cases of chaotic dynamics seems to limit seriously its applications. Despite this one can 
try to understand with it the scar phenomena noticed by Heller [14] and investigated by the 
latter author and others [15, 19, 20, 21]. Although in examples discussed below we do not 
construct closed skeletons but yet supposing their existence in these cases allows us at least 
half-qualitatively to understand the scar phenomena and to make some predictions for their 
existence in some cases of billiards and their absence in others. 

Consider for example the Bunimovich stadium billiards of Fig. 10 and the isolated unstable 
horizontal periodic orbit linking the top points A and B of the stadium. It is known [14] that 
the there is a mode for which its wave function takes significantly larger values around the 
orbit than far away of it, i.e. the orbit signals its existence by such a "scar". This and 
similar scar phenomena have not got their full description although there were many efforts 
and approaches to do it [15, 19, 20, 21]. 

Of course it is not easy not only to construct a skeleton carrying such a mode but even 
to prove its existence. Nevertheless assuming the latter we can expect that such a skeleton 
will be symmetric horizontally and vertically as well as its bundles will contain the horizontal 
periodic orbit. 

Consider one of such bundles emerging from a vicinity of the top point A on Fig. 10. Even 
if its rays are divergent with respect to the periodic orbit the rays which are very close to it 
are also almost parallel to it, i.e. after the reflection by the opposite semicircle of the billiards 
they are transformed into convergent bundle focused close to the focal point of the reflecting 
semicircle. It means that central parts of almost all such bundles of the skeleton containing 
the horizontal periodic orbit have to be convergent and have to pass close to focal points of 
the semicircles. 

Assume the radii of the semicircles to be equal to 1 while according to Fig. 10 the flat 
parts of the stadium have the length a each. It is clear that each such a convergent bundle 



?/, A) = Asm{Xpix) sm{Xp2y) 




\pia = rmT, Xp2b = mr, m,n = 1,2,... 



(128) 



36 



will generate after the reflection a new bundle the central part of which is again convergent. 
On this new bundle a new x-fector defined on it will however be weakened according to (40) 

— - 1 

by the factor 2 ^ _|_ 3^-5. Therefore if the starting bundle has the factor xo 

as its zeroth order semiclassical approximation then after n subsequent reflections it will be 
weakened by the factor (2a + 3)~2". However the central parts of all the reflected bundles 
will remain close to the periodic orbit so that SWF's defined on them will interfere in an 
infinite number of them close to the orbit. 

Easy calculations of the contributions coming from all such bundles lead us to the following 
form of the regular scarring part of the semiclassical wave function on the horizontal periodic 
orbit of Fig. 10 in the JWKB approximation: 

^JWKB( ^-^ ^ 
^ scar y-^i 

piXp{a+2) ^ J 

1 _ g2,2ap(a+2) ((2- + 1)-^ e-^-(— - (2(a - x) + l)"^ ,e+^A.(a-.+l)) ^^(1^ 

-l<x<a + l, ar^-^,a + ^, g = (2a + 3)-^ (129) 

where arg(2x + 1) = vr for .t < — ^ and the point x = — ^ is avoided clockwise and above it 
when X moves to the values larger than — | and the point A has been chosen as the initial 
one for the s-variable so that s y — 1 in a vicinity of this point. 

Certainly the above contribution is not a unique one (^'^^^^(x, 1) does not vanish for 
X = — l,a + 1, see our comment below). There are still infinitely many bundles contributing 
to the SWF coming from the parts of the bundles discussed above scattered away of their 
central parts. However it seems reasonable to assume that just such bundles scattered upon 
the whole billiards interfere everywhere chaotically, i.e. typically as in the cases without scars 
and satisfying Berry's conjecture [26]. Some order in such interfering can be organized as we 
can try to argue by the presence of periodic orbits supported by the focusing properties of 
the billiards boundary in vicinities of top points, i.e. the points where periodic orbits touch 
the billiards boundary. 

The last condition, i.e. the focusing properties of the billiards boundary at the top points 
seems to be essential for the scar phenomenon to appear. One can be convinced of its 
necessity by considering an " anti-Bunimovich" stadium, i.e. the stadium which semicircular 
parts instead of being concave are convex for the billiards, see Fig. 12. The horizontal periodic 
orbit still exists but a skeleton with the properties discussed earlier certainly does not, i.e. 
the bundles containing this periodic orbit can be only scattered so that the only rays of these 
bundles which come back close to the periodic orbit is the orbit itself. All the remaining rays 
even close to the periodic orbits are scattered away passing by the virtual focus points of the 
convex semicircles. 

The results of Barnett's calculations demonstrated by Sarnak [24, 25] seem also to confirm 
this conclusion. 

Let us note however that the above analysis is not sufficient to estimate the energy cor- 
responding to the mode described above even its lowest JWKB approximation investigating 
only a vicinity of the horizontal periodic orbit . This is because after every reflection the 
corresponding half of the horizontal periodic orbit belongs to a new bundle, i.e. this orbit 
is in fact not closed on the skeleton considered and moreover it does not end in the initial 
bundle (this is why it is not periodic inside the skeleton). Because of that this classically 
periodic orbit cannot be used to write on it "the last quantization condition" of the form 
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Figure 12: The anti-Bunimovich billiards with the bouncing mode skeleton deprived of the 
scar generating bundles associated with the horizontal periodic orbit 



(49) . In fact such a single periodic orbit possibility of determining the energy would question 
the role of the eigenfunction boundary conditions in forming its corresponding eigenvalue. 

9 Summary and conclusions 

We have shown in this paper that it is possible to formulate the semiclassical description 
of the quantum billiards eigenvalue problems in the semiclassical wave function language. 
We have used in principle the approach of Maslov and Fedoriuk [6] modifying it however 
by the way of moving through caustics. The corresponding procedure has been described in 
App.B. where the idea of continuing semiclassical wave functions by coustic singularities on 
the complex t-plane was developed. It was shown there a close relation between a signature 
of a SWF and a path along which the SWF has to avoid a caustic singularity on the t-plane. 

We have also modified the Maslov and Fedoriuk procedure by constructing respective 
Lagrange manifolds not as smooth, deprived of boundaries tori-like surfaces but rather in the 
form of skeletons. The skeleton idea has appeared to be sufficiently flexible to gather quantum 
systems at least in principle uniformly independently of the kind of their semiclassiclal limits, 
i.e. whether these limits are integrable or chaotic. Nevertheless in the last cases constructions 
of the corresponding skeletons seem to be rather difficult except these obvious chaotic billiards 
boundary configurations which have been considered in this paper. 

It is worth to stress also at this summary that the forms of skeletons as a set of ray bundles 
satisfying the geometrical optics law of the mirror-like reflection from the billiards boundary 
was not a matter of a free choice but a necessity of satisfying the boundary conditions by the 
SWF's as it was shown in App.A. 

No less important was establishing that in the semiclassical calculations the zeroth order 
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Figure 13: The example of the four bundles skeleton in the right angle triangular billiards on 
which none semiclassical eigenfunction can be defined 



terms of SWF's are classical integrals of motion. This has allowed us to close the correspond- 
ing calculations. 

Using known examples of the classically integrable billiards such as the circular and 
rectangular ones and the broken rectangular billiards as some variants of the rectangular 
ones we have demonstrated the effectiveness of the skeleton method in the quantization of 
the systems mentioned. In particular the semiclassical calculations of energy levels performed 
in sec. 5 have developed a new algorithm of finding these levels to any order in Planck constant 
for the case of the circular billiards (i.e. the cylindrical infinite well). 

Applications of the method to other classically integrable billiards systems also in higher 
than two dimensions are in preparations. 

We have shown also how the skeleton method allows us to describe almost trivially the 
bouncing ball modes in the Bunimovich billiards and similar ones. 

We have also discussed a possible explanation of the scar phenomena by the skeleton idea 
showing a close relation of these phenomenta to focusing properties of billiards boundary in 
vicinity of the top points of the periodic orbits carrying scars. A lack of focussing properties 
by the billiards boundary in vicinity of the top points should exclude the existence of scars 
around the corresponding periodic orbits. 

On the other hand it should be stressed that constructing a desired skeleton does not 
guarantee a construction on it a SWF satisfying desired necessary boundary and other condi- 
tions. A simple example of such a situation is a rectangular triangel billiards shown on Fig. 13. 
There is a skeleton collected of four bundles as shown in the figure. It is clear however that 
none SWF vanishing on the boundary of the billiards can be built on this skeleton. 
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Appendix A 

In this appendix we are going to show, that the geometrical optics rule of reflections of 

rays off the billiards boundary is a consequence of demands of vanishing on the boundary of 
the linear combination (35) accompanied by the conditions (36) and (37). Namely consider 
the following superposition of SWF's: 

^l'{x,y;u,l;\) = ^l]-^{di, sv,u,l; X) + ^l^^{d2, S2;u,l; X) = 
J'l (di, si;u, /)e''^iH<ii+/;' cosa,,i(.';«,0<i.')4i^(di, si-u, I; A) + 

J~l {d2, S2;u, l)e''''K'^'+fo' ^°^'^'=-2(^>'')'^^');^^22(d2, S2; u, I; X) (130) 

with 

r = [x,y] = rk,i{di,si;u,l) = ro(si) + di{si;u,l) = 

rfe,2(o!2, •S2; u,, I) = ro(.S2) + d2(.S2; u, I) 
rk,i{d, s; u, I) € Bkiu, I), Vk,2{d, s; u, I) € B'j^{u, I), Bk{u, I) / B',^{u, I) (131) 

i.e. the SWF's '^^\{di,si;u, I; X) and '^^\{d2, S2\ u, I; X) are defined respectively on the bun- 
dles Bk{u, I) and B'f,{u, I) with Dk{u, I) H D'f^{u, I) interfering in the crossing point [x, y] 
of two rays r^. i(di, si; n, /) and 2(t^2) ■52) belonging to the respective bundles. 
Therefore the condition for \l'^*(x,y, A) to vanish on Ai^{u,l) is: 

Jk} (0' ^y'''" cos akAs';u,l)ds'^a^^^^. ^ 
Jk,l (0' ^y''"'^ ^° cosak,2{s';u,l)ds' ^.2^^^. A) = 

ro{s) e Ak{u,l) (132) 
Because of the fc-dependence the last relation can be satisfied if and only if: 

rs rs 

a\ \ COS ak,i{s' ; u, l)ds' = a2 / cos ak,2{s' ', u, l)ds' , ro{s) E Ak{u,l) (133) 
Jo ' Jo ' 

It is easy to see however that there are only two solutions of the last condition: 

ak,i{s;u,l) = ak,2 {s;u,l) for ai = 02 
ak,i{s;u,l) = Tr - ak,2{s;u,l) for ai = -a2 

ro{s) e Ak{u,l) (134) 

The first solutions are however uninteresting identifying the bundles in a given segment 
and consequently leading to the solutions vanishing identically on Ak{u,l). 

Putting ak,i{s;u,l) = a{s;u,l) and ai = —a we get from the second solution and from 
(132): 

Xk^{s; u, I; A) = -Xk,i{s; u, I; A) = Xk{s; u, I; A) 

ro{s) e Ak{u,l) (135) 



40 



so that the combination (35) becomes: 

^I'ix, y; u, I; A) = *^_i(di, si; u, I; A) - ^^^"2(^2, S2; u, I; A) = 

J-| id2,S2;u, l)e-'<''+Io' -^"(^'^'^'')'^^');^-(d2, S2; u, I; A) 

ro{s) e Ak{u,l) (136) 

where XkO'id, s; u, I; A), <t = it, are given by (31) with Xfc(O) s'l u, I; A) = Xk{s', u, I; A). 

The last result shows that '^'j^^{x,y;u,l; X) vanishing on Ai^{u,l) has to be represented 
semiclassically by a combination of at least two SWF's of opposite signatures and such that 
if '^^i{d,s;u,l;X) is defined on the bundle Bk{u,l) then the second SWF ^'^2{d,s;u,l; X) 
has to be defined on the bundle B^{u, I). 



Appendix B 

In this appendix we explore the well known properties of the circular billiards to establish 
the way of avoiding the singular caustic points on the t-plane as well as to confirm the way 
given by (38) by which two SWT's have been matched. 

The circular billiards is of course the well known case of the two dimensional infinitely deep 
cylindrical potential well which energy spectrum is easily obtained by solving the stationary 
Schrodingcr equation (SSE) for this case by the variable separation method performed in the 
cylindrical coordinates. In the latter coordinates the radial part of the SSE is the following: 

^"(r) + ^^'(r) + X^(^E- ^ j ^(r) = (137) 

where the separation constant m = 0, ±1,..., is the angular momentum quantum number, 
A = Ti~^ and we have put 2M = 1 where M is the billiard ball mass. 
After the substitution *(r) = we get: 

r2 



<(r) + X^E- V(r) = (138) 




AV 

Assuming the unit radius of the billiards we get the solutions for (137) in the following 
forms [27]: 

'^±{r) = g~^e^'^-^i'^^*'x±(^), r+<r<l 
q = p, = E-^ = --^ir-r+){r-r-), r± = ±j-j=, m/0 (139) 

which correspond to the fundamental solutions ip+{r) and tp-{r) defined in the respective 
sectors 1 i I of the Stokes graph of Fig. 13 with x±{r) normed at the corresponding sector 
infinities by x+(ooi) = x_(ooj) = 1. The r-plane on this figure has a cut from r+ to r_ 
on which changes its sign, while gets the factor ±i depending on the cut crossing 
directions {+i for crossing it upwards, —i - for downwards). 



41 




1 



Figure 14: The Stokes graph for the cyhndrical bihiards 



The Stokes graphs formahsm corresponding to (139) easily provides us with the energy 
spectrum of the bilhards. For this we take the wave function ^{r) corresponding to this case 
as a linear combination 

^(r) = ijj+{r) + cijj-ir) (140) 

given on the segment AC where c is defined by the condition V'(l) = 0- 
It is easy to see that then: 

c=-¥^, (141) 

Next the combination (140) has to be continued to the segment OC vanishing on it 
when r — )• 0. This means however that up to a constant it should be identified with the 

fundamental solution ipo{r) = q~ie •^'■+ Xo('')) taken for example in this form above 
the cut OC attached to the sector on Fig. 13, i.e. for r + ie, < r < r+. 

According to the well known rules the latter can be expressed as the following linear 
combination of ip±{r) [27]: 

~ iX qhdr' ~ -i\ qhdr' ~ 

< r < r+ (142) 
where Xo^j = limj._i.ooj Xoi^) with ooj being the infinity point in the sector j, j = 1, 1. 
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It then follows from (142) that: 



rl 1 

tjj{r) = e ■^'■+ V'o(r) (143) 

Xo->i 



and 



c = I— e + (144) 

Xo-)-i 

The equations (141) and (144) together are just the quantization conditions for the energy 
with the spectrum of the latter being degenerate with respect to the sign of m. 

To make further a correspondence with sec. 5 we shall use rather the solutions ip±{r) = 
X=f(1)V'±('")) X±(^) = XTi'^)x±iT) and il){r) = X-{^)'^{f) in the remaining considerations. 
Therefore instead of (140) we get simply: 

V'(r) = ip+{r) - V'-(r) (145) 

so that and instead of (144) we have: 

rl 1 

Yi_s.n -2iX q'^dr' 

1 = -iii±ie ■''■+ (146) 

where Xi^o = X-(l)Xi-)-o and Xi^o = X+(l)Xi->o and where we have taken into account 
that xo^j = Xj^o, j = 1, i- 

The results (141)-(146) obtained above are exact. 

Consider now the semiclassical limits of (141)-(146), i.e. when A +oo. It should 
be stressed that the fundamental solution formalism allows us to do such a passage and in 
every of the above formula this limit is well defined. In particular all exponentially small 

contributions have to be neglected. In this way if r is in sector such an exponentially small 
is the solution ^o{r) so that taking the semiclassical limit in (142) we get: 

1 ^ 1 ^ 

= -zXi!.oe^'''^^+ (0 + X?i.oe"'^^+ '^''"^'rJir), < r < r+ 

J 1 

r± (r) = Qa? e^'^ /i" «-*'x5? (r) (147) 

where "as" denotes the semiclassical forms of the relevant quantities and in particular 
X±{r), xff^o(~ Xo^i) and the asymptotic form of Xi^o(~ Xol>i) t)e found in any of 
the references [27] while qas = E — with E = Eq + ^ .... 

In fact the linear relation (147) is equivalent to the following identity valid in the semi- 
classical limit: 

X+{r)xlU^ = xlUoX-{r), < r < r+ (148) 

if the point r, r G OC, is approached continuing x+{i^) anticlockwise and x"*(r) - clockwise 
around the point C of Fig. 13. 

Taking also the semiclassical limit of the exact quantization condition (146) we get: 

X?f^n -2iA q^dr' , , 

l = -^^e (149) 
Xi->o 
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Prom both (147) and (149) it follows that: 



iPl'ir) = 'iljl'{r), 0<r<r+ (150) 

Conversely if (150) is valid then from the identity (147) we get the semiclassical form 
(149) of the quantization condition. 

Therefore both the conditions (146) and (150) are equivalent to each other as the semi- 
classical quantization conditions in the case considered. 

We conclude therefore that the equation (150) substitutes the quantization condition (149) 
in the semiclassical limit, i.e. in the quantized semiclassical limit both ijj'^{r) and il}"f{r) has 
to coincide in the sector 0. But this coincidence can be analytically continued to any (non 
singular) point of the r-plane. In particular it can be continued back to the classically allowed 
segment AC. These can be done in two ways, by avoiding the singular point C clockwise or 
anticlockwise. 

Before making the respective continuations let us note that the relation (142) cannot be 
continued to the segment AC of the Stokes graph of Fig. 13 not loosing its connection with its 
asymptotic form (147) while this asymptotic form itself can loosing however its connection 
to the original equation (142). This new form of an equation substituting the original one 
(142) can be obtained by the Borel resummation procedure. 

Let us start therefore with the relation (147) and (150) to continue them to AC around 
the point C clockwise. We get: 

1 J 1 

(151) 

and 

V'+ (r) = V- ''°"*(r), r+<r<\ (152) 

where il;°f''^"'*(^r) denotes ip°f{r) continued analytically clockwise around the point C from 
the segment AC to the same segment. 

If we choose the opposite path of continuation, i.e. anticlockwise one, then continuing the 
relation (150) we get: 

1 ^ J 1 

(153) 

and 

^«^(r) = V'"' ''°"*(r), r+ < r < 1 (154) 

where ip°^''^°"'^(r) denotes ^^^^(r) continued analytically anticlockwise around the point C from 
the segment AC to the same segment. 

However, contrary to OC the segment AC lies now totally in the common part of the 
domains of Borel summability of iIj± (r) so that we have: 

(V'± (r))^^ = V±(r-)> r+<r<l (155) 
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where the superscript BS means the Borel resummation of the relevant quantity. 
But according to (152) i (154) we have instead: 



r+ < r < 1 



(156) 



while by the Borel summing (151) and (153) we obtain respectively: 



(I as.cont / \\ 
■4>- in ] 



BS 



(I as.cont/ \\ 
^+ [r) ) 



BS 



I e 

Xi-i>o 



V'-(r) 



< r < r+ 



(157) 



If we now apply the Borel summed quantization conditions (156) to the respective equa- 
tions in (157) we immediately recover the exact quantization condition (146). 
Let us summarize the above results. 

• In the semiclassical limit the quantization condition for the circular billiards is the 
appropriate coincidence of the asymptotic expansions of two fundamental solutions in 
the classically unallowed region (the segment OC). This coincidence is given by (150). 

• In the classically allowed region the semiclassical quantization condition can be obtained 

by continuing appropriately (clockwise or anticlockwise) one of the semiclassical solu- 
tions ^p'jf{r) from the segment AC and back to it around the turning point C and next 
by identifying the continued solution with the second one. The respective identifications 
are given by (152) and (154). 

• The semiclassical quantization conditions formulated in the classically allowed as well 
as in unallowed regions can be Borel summable leading us to the exact ones; 

• The fundamental semiclassical solutions continued to the classically allowed region (the 
segment AC) by the unallowed one (the segment OC) still represent semiclassical ex- 
pansions of some fundamental solutions which by Borel resummation can be identified 
appropriately (see the formulae (157)). 

• The semiclassical quantization condition formulated by (152) and (154) in the classically 
allowed region (the segment AC) can be Borel summable providing us with the proper 
exact quantization conditions. 

B.l The fundamental solutions semiclassical analysis rewritten on the clcis- 
sical trajectories 

Let us now rewrite the above results in the formalism of the SWF's defined on the classical 
trajectories developed in sec. 5. 

To this goal let us assume the semiclassical expansion (12) for the energy E so that the 
classical momentum p = \/2£o) and associate with the solution (140) two bundle skeletons 
B and B^. 
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Both the skeletons consist of single ray bundles only. The skeleton B consists of the 
ray bundle .6(0, 27r) the rays of which run with the definite positive angular momentum 
m, m = 1,2, and with momenta p shown in Fig. 2 while the angular momentum of the 
skeleton B"^ is equal to —m and rays of its unique bundle B^{0, 27r) have momenta p"^ on 
Fig. 2. For both the momenta we have |p| = jp^^l = p. 

Any point of the ring (r, r_)_ < r < 1, < $ < 27r, is crossed by two rays coming 
from different bundles starting from the points si and S2 of the billiard boundary and having 
the respective angular momenta equal to pcosa = m/X and pcos(7r — a) = —m/X with 
cos a = r-i- (see Fig. 2) while their corresponding radial momenta are the same and equal to 

/ 771 cos^ a 

-y^o f2 — ■ 



r{t) 



^ = -q^ being negative {t < tc, tc 



sin a/p). 

Starting at the moment t = from the boundary these rays meet each other at the point 
(r, $) in the moment t. 

Each of the bundles defines a transformation of coordinate (r, $) — )■ {d,Si), i = 1,2, by 
r = r{d) = yjcos^ a + (sin a — d)^, d = pt and $ = si + a — 7 



correspondingly to the bundle, see Fig.2. 
The total original solution 



1 
f 2 



S2 — a + 7 < si + a 



(158) 



corresponding to (145) can be now rewritten by the new coordinates noticing that q^dr' 
- jlprdr' = - j^pldt' = - J^p'^dt' + /o* P^dt' = -pH+p-" cos2 a and that 



* dt' 



r^(t') cos a 



arctan 



p{t' - Q 



cos a 



cos a 



(a -7) 



(159) 



for t <tr so that: 



x±it,x) 



^±{r{t),X) = go *e^^^P'*±*l-l("-^)x±(t,A) 



X±{rit),X) 



— exp 



Ti^^n+iA-"-^ r ^ -dt' 



1 + 

Q_ 



1 + 



X±{r{t),X) 



Qo 



-n-l 



(160) 



and then we get: 



*(r,$) = 



V>+(r(^))-V;_(r(t)) 

T^^ ' 



(r2go)"^e-^^P'*+'l'"l("-^)+^"**X+(i, A) - (r2go)-3e+^^P'*-^l"'l("-T)+^"^*X-(i, A) (161) 
where {r^qo)~^ = {r^{t)qo{t))-^ = {Eo{r'^{t) - rl))'^ = p-^\t - tcl"^ . 
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Prom (161) we get further: 



{r^qo)-h-'^P''*+'"''^X+it, A) - (r2go)"^e^^^''*+^"*^2^_(t, A), m < 

f i^"(t,S2) - m>0 
1 - i*+(t,S2), m<0 



(162) 



where ^"^{tjSj), j = 1,2, have the form (26) but are exact. We have also assumed that 
{t — tc)^2 = —i\t — tc\~^ for t < tc (see below). 

The last result shows that, depending on m, the bundles on which the fundamental 
solutions '4>±{r) are defined are chosen by the latter accordingly to their signatures. 

Let us now consider the analytical continuation of the semiclassical limit ^'"^(r, ^>) of the 
solution ^{r, $) from the point A to its point B as it is defined by the rule (152) or back, i.e. 
from B to A correspondingly to the rule (154), see Fig. 2. 

Since t = tc — | -\/(r — r+)(r + r+) then, if r rounds the turning point C by 2it clockwise 
starting from the segment AC (see Fig. 13), t rounds tc on the t-plane also clockwise by vr. 
But t rounds tc by tt anticlockwise if r does the previous motion in the opposite direction. 
In both the cases we are found at the point B with t > tc- However in the first case 
{t — tc)~^ = \i ~ tc\~^ while in the second {t — tc)~'^ = —\t — tc|~^. 

^"■^{r, $) continued clockwise on the r-Riemann surface is given by: 

KLkir, $) = ^+^^'^';-'^^""^'^ e-^ = 4 (l - ^e-'''^^' '''''' ^) (0^^"^* 

r2 r2 V Xi-*-o/ 

r+ < r < 1 (163) 



where ip°f''^*(^r) has been substituted by its form which follows from (153). 
Continued however anticlockwise it is given by: 

r+ < r < 1 (164) 

where ijj"^''^'°'^^{r) has been substituted by its form given by (151). 

Note that in both the above continuations we have not taken into account the quantization 
conditions (152) and (154). 

However if the energy is to be quantized we have to have: 

^cLkir, ^) = K^ticlockir, ^ (165) 

Performing now in (163) and (164) the calculations similar to the ones above but with 
the variables t, t > tc, we get from (162) and (163), see Fig. 7: 

VI/- . ^.^f -i^-^^'"''\t,si)+W-{2tc-t,S2), m>0 ^ 
clockK'^ ) \ -i^+''=°'^*(i,s2) + i*-(2tc-i,si)), m<0 
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P~^\t - t^|-5e-^V(2te-t)-iAps2Cos(7r-a)^as(2t^ -t,X), m > 



. - i^|-2e-*^P'(2t,-t)-iApsicosa^as(2tc - t, A), m < 

2tc>t> tc, S2 = si± 2a (166) 

where s^), f > t^, are the results of the continuations of *"*"(*, sj), t < tc, j = 1, 2, 

by the caustic. 

For Knticlocki''^^'^) we get instead: 



-i^'+(2tc - t,si) +i^'"'™"*(t,S2), m>0 _ 
-i^+(2tc -t,S2) +i*"'''°"*(i,si)), m<0 

-p-^\t - t^|-|e*^P'(2tc-t)+iApsicosQ^as^2tc - t,X) + 

-ip-^\t - f^|-5^e-*^P'*-^^P"2COs(7r-a)^as(^^^)^ ^ > g 

_ ^^|-2eiAp2(2t,-t)+iAps2Cos(7r-a)^as|-2^^ _ t,\) + 

-ip-^\t - tc\-^^e-'^P^*-'^P'^''''^'^X- {t, X), m<0 

2tc>t> tc, S2 = si± 2a (167) 



where s^), f > t^, are the results of the continuations of '^"{t, sj), t < tc, j = 1, 2, 

by the caustic. 

From the above considerations it follows clearly that ^"""(t, s) continued through the caus- 
tic has to avoid the focusing point at t = tc from above in the t-plane, moving clockwise while 
'^~{t,s) - from below, moving anticlockwise. By such a continuation *+(t, s) acquires the 
factor i while ^^{t, s) - the factor — i. 

Let us summarize the above results. 

1. The SWF 'I'"*'(.x,y) for the circular billiards originally defined closely to the boundary 
of the billiards and vanishing on the boundary can be written as the combinations (162) 
of the two SWF ^^{t, s) of the opposite signatures having the forms given by (26) and 
defined on the two skeletons B and B"^. 

2. It is the SWF ^'^(t, s) which is defined on the skeleton B which momenta being tan- 
gential to the caustic of the bundle give the anticlockwise orientation of the caustic for 
m > but on the skeleton B"^ for m < 0. The SWF '^~(t,s) is then defined on the 
second skeleton in the respective cases. 

3. ^""^{x, y) can be continued across the caustic avoiding the real singular point at t = tc 
on the t-plane from above for ^~^[t,s) and from below - for ^~(t,s) and getting then 
the forms (166) or (167). 

4. If \I'"*(x, y) is to be quantized then ^«^.cont|-^^ j^^g vanish identically independently 
of the way (clockwise or anticlokcwise) it was continued. In such a case we have to have: 

ii,±,'^t(t,s) = '^^{2tc-t,s±2a), t>tc (168) 
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leading to for t = 2tc- 



exp 



J2 En+iX-"" / -dt' 



n>0 



^iXp^2tc—2iXap cos a+i^-ir _ (169) 



where we have taken into account that r{2tc — t) = r{t). 
5. In the hmit A ^ oo the condition (169) decays into the following ones: 

^iXpd{s,a)—2iXap cos a+i^TT j 



Xi!.o(^,A) 
xr-.oiE,X) 



cxp 



i J2 ^n+lA-" / -dt' 

n>0 1 + A V 



d{s, a) = 2ptc = 2 sin a 



(170) 



Of course the condition (169) coincides with (144). 

6. The two degenerate SWF's \I'"*'^(x, y) in each point inside the ring cos a < r < 1 can 
be given by: 



^r«^'+(a;, y) = -i^+(t, Si) + i^+''"'''\2tc -t,S2- 2a) = 
i^-{t, S2) - i^-''''"'\2tc -t,si + 2a), m > 
^r«^'-(a;, y) = -z^'+(t, S2) + i^+''^\2tc -t,si- 2a) = 
_i^-,cont^2tc -t,S2 + 2a) + i^-{t, si), m < 

t < tc 



(171) 



(172) 



Appendix C 

It is shown in this appendix that 6k{u,l) from the formula (33) and 5(uj,lj) from the 
formula (39) are s-independent. To this end consider Fig.l on which a mapping hk{s;u,l) of 
the arc Ak{u, I) into an arc Ak'{u', I') is defined by the the bundle Bk{u, I). According to this 
mapping we have: 



ro{hk{s; u, I)) = ro(s) + D(s; u, I) 
ro(s) G Ak{u, I), ro{hk{s; u, I)) G Ak'{u', I') 



(173) 



where D(s; u, I) is a vector linking the point ro(s) with the point ro(/ife(s; u, I)) of the billiards 
boundary. 

Differentiating the equation (173) with respect to s wc get: 

dhk{s;u,l) djkis;u,l) . 

cos p{hk{s;u, I)) = cosp{s) - D{s;u,l) sm^k{s;u,l) + 



ds 



ds 

cos7fc(s;ti,/) 



dD{s; u, I) 
ds 



dhk{s;u,l) . • N , r>r ,^d'yk{s;u,l) 

■smp{hk{s;u,l)) = sm^(s) + D{s;u,l) cos7fe(s; u, <) + 



ds 



ds 

sin7fc(s;ti,/) 



dD{s; u, I) 
ds 



(174) 
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where D{s; u, I) is the length of D(s; u, I). 

Making further the proper linear combinations of the last equations we have finally: 

dhk{s;u,l) , , dD{s;u,l) 
cos ak'{hk{s;u,l);u ,1 ) = cos ak{s-u,l) H 

OS OS 

smak'[hk[s;u,l);u ,1 ) = smak[s;u,l) - D[s;u,l) — (175) 

where we have taken into account the following relation between the angles involved: 

ak'{hk{s; u, l)]u', I') + Ukis; u, I) = (3{hk{s; u, I)) - /3(s) + 27r 

Jk{s;u,l) = P{s) + ak{s;u,l) (176) 

which follows from Fig.l. 

The independence of s of Sk{u,l) and S{uj,lj) follows now easily from the first of the 
relations (175). 

Appendix D 

The billiard Laplacean expressed by the variables d, s has the following form: 

/ cos^ a .\ d 19^ 2 cos a 5^ 
A{d,s) = y—^ + + J^d^ ~ J2 dsdd ^ 

fcosa [. , y . \ 1 /7" , . \ j'\ d 

[-J3- (^(a +7 ) cosa - — smaj + cosa + a sma j - -j — - 

-^(^(a+7)cosa--smaj+--^j- 

J = sina — 7'd (177) 

while for the corresponding operator J2 • A • J^2 we get: 

1 1 3 1 

J5 . A • J-5 = j-^(yjf - j-^yj . V - -J-^ A J + A 
A ^ ' 2 
/2 2 ' 1 

iyjf = '^{S - 2d sin a + l) + cos a - d-y") cos a + j^{a' cos a - d'y"f 

V J ■ V = — ^-^2 {d^ — 2(i sin a + 1 j ^ ~ '^'^"^ "^^^ ^ ^d'^ 

f 1 , , , 7' cos a \ d 

(^-ia'cosa-di') + :L__^_ 

15%. , 27"cosa 
^■^= j2^(sina-7ti) + — j2 

-j^ [{a cosa — ^7") cos a — (/3' sino — d'y''^){sm.a — j'd)^ — ^ — (a' cos a — dj") (178) 

For the circular billiards the corresponding Laplacean takes the forms: 

9^ 1 d"^ 2 cosa d'^ Id cosa d 



Aid, .) = (d2 _ 2dsin a + 1) + -1-^ 



J2 V J dd'^ J2 J2 dsdd J dd J3 9s 

J = sinQ;-d (179) 
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while the corresponding operator J2 ■ A ■ J 2 the form: 

tIatI/.n 3/,9^,. \ l/,9„,. \ d cos a d 

J2 . A ■ J'-2{d, s) = ^ [d^ - 2dsma + l) + [d^ _ 2dsina + ij — - + 

^+A(d,s)(180) 

For the rectangular billiards the corresponding Laplacean has the following two similar 
forms depending on the rectangular sides: 



1 / 9^ d"^ \ 2 cos a d' 



,2 



^^^'^^ sin^a ydd? ds^ J sin^a 5s5d 

_ 1 / 9^ ^\ 2 sing d'^ 
cos^a ISc/^ ds'^ J cos'^ a dsdd 

and of course J2 • A • J~2 = A for this case of the billiards. 
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